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frt The present review analyzes and compares various types of dissipationless spin transport: (1) 

^^ Superfluid transport, when the spin-current state is a metastable state (a local but not the 

^^^ absolute minimum in the parameter space). (2) Ballistic spin transport, when spin is trans- 

ported without losses simply because sources of dissipation are very weak. (3) Equilibrium 
C^ spin currents, i.e., genuine persistent currents. (4) Spin currents in the spin Hall effect. Since 

superfluidity is frequently connected with Bose condensation, recent debates about magnon 
Bose condensation are also reviewed. For any type of spin currents simplest models were 
chosen for discussion in order to concentrate on concepts rather than details of numerous 
models. The various hurdles on the way of using the concept of spin current (absence of the 
spin-conservation law, ambiguity of spin current definition, etc.) were analyzed. The final 
conclusion is that the spin-current concept can be developed in a fully consistent manner, 
and is a useful language for description of various phenomena in spin dynamics. 
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Introduction 



^ The problem of spin transport occupies minds of condensed matter physicists for 

■^ decades. A simple example of spin transport is spin diffusion, which is a process 

CN accompanied with dissipation. Conceptually more complicated is "dissipationless" 

_J spin transport, which was also discussed long time but was in the past and remains 

now to be a matter of controversy. The main source of controversy is that spin is not 
Q a conserved quantity. This leads to many complications and ambiguities in defining 

00 such concepts as spin flow, current, or transport. Sometimes these complications 

^^ are purely semantic. However, this does not make them simpler for discussion. 

^ "Semantic traps" very often are a serious obstacle for understanding physics and 

for deriving proper conclusions concerning observation and practical application of 
the phenomenon. The best strategy in these cases is to focus not on names but 
^ on concepts hidden under these names. Only after this one may "take sides" in 

semantic disputes not forgetting, however, that choosing names is to considerable 
extent a matter of convention and taste. 

During long history of studying the problem of "dissipationless" spin transport 
one can notice three periods, when studies in this field were especially intensive. 
The first period started from theoretical suggestions on possible "superfluidity of 
electron-hole pairs" [IJ, which were later extended on possible spin superfluidity 
[21 13] . At the same period the concept of spin superfluidity was exploited j^^ for 
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interpretation of experiments demonstrating unusually fast spin relaxation in '^He- 
A [7J. The second period was marked by intensive theoretical and experimental 
work on spin superfluidity in ^He-B starting from interpretation of experiments 
on the so-called Homogeneously Precessing Domain (HPD) j8] in terms of spin 
supercurrents [9\. Finally in these days (the third period) we observe a growing 
interest to dissipationless spin currents in connection with work on spintronics 
|10j . The final goal of spintronics is to create devices based on spin manipulation, 
and transport of spin with minimal losses is crucial for this goal. Now one can find 
reviews summarizing the investigations done during the first [llj and the second 
|12^ [T3] periods of works on dissipationless spin transport. On the other hand, 
the work on spin transport in spintronics is a developing story, and probably it is 
still premature to write summarizing reviews. Nevertheless, some reviews mostly 
addressing the spin Hall effect have already appeared [HI [E] . It looks also useful to 
have a glance on the current status of the field from a broader viewpoint and to find 
bridges between current investigations and those done in the "last millennium". 
The present review aims at this goal. The intention is to discuss mostly concepts 
without unnecessary deepening in details, and simplest models were chosen for this. 

The term "superfiuidity" is used in the literature to cover a broad range of phe- 
nomena, which have been observed in superfluid "^He and ^He, Bose-Einstein con- 
densates of cold atoms, and, in the broader sense of this term, in superconductors. 
In the present review superfluidity means only a possibility to transport a physi- 
cal quantity (mass, charge, spin, ...) without dissipation. Exactly this phenomenon 
gave a rise to the terms "superconductivity" and "superfluidity" , discovered nearly 
100 years and 70 years ago respectively. It is worthwhile to stress that one should 
not understand the adjective "dissipationless" too literally. In reality we deal with 
an essential suppression of dissipation due to the presence of energetic barriers of 
the topological origin. How essential suppression could be, is a matter of a special 
analysis. In the present review we restrict discussion with the question whether 
activation barriers, which suppresses dissipation, can appear. 

But superfluidity is not the only reason for suppression of dissipation in the 
transport process, and it is important to understand the difference between various 
types of dissipationless transport. In the present review we shall discuss four types 
of them: 

• Superfluid transport: The spin-current state is a metastable state (a local but 
not the absolute minimum in the parameter space). 

• Ballistic transport. Here spin is transported without losses simply because 
sources of dissipation are very weak. 

• Equilibrium currents. Sometimes symmetry allows currents even at the equilib- 
rium. A superconductor in a magnetic field is a simple example. Equilibrium spin 
currents are also possible, though there is a dispute on whether they have some- 
thing to do with spin transport. Equilibrium spin currents are genuine persistent 
currents, since no dissipation is possible at the equilibrium by definition. 

• Spin currents in the spin Hall effect. These currents are also called dissipationless 
since they are normal to the driving force (electric field) and therefore do not 
produce any work. However, in the spin Hall effect there is dissipation connected 
with a longitudinal charge current through a conducting medium. On the other 
hand, it was recently revealed that the spin Hall effect is possible also in insulators 
where a charge current is absent. Then spin currents are not accompanied by 
any dissipation becoming similar to equilibrium spin currents. 

The second type (ballistic) looks mostly trivial: dissipation is absent because 
sources of dissipation are absent. Still it is worth of short discussion since some- 



March 3, 2010 17:4 Advances in Physics SpinRev 

Advances in Physics 3 

times they confuse ballistic transport with superfluid transport (an example of it 



is discussed in section 7.2). The superfluid transport does not require the absence 
of dissipation mechanisms. One may expect that in an ideally clean metal at zero 
temperature resistance would be absent. But this would not be superconductivity. 
Superconductivity is the absence of resistance in a dirty metal at T > 0. 

The first two types of spin currents are discussed in Part I of the review, which is 
devoted to magnetically ordered systems. The third and the fourth types are dis- 
cussed mostly in Part II, which addresses time-reversal-invariant systems without 
magnetic order, though in magnetically ordered media equilibrium spin currents 
are also possible (section^. Since from the very beginning of the theory of super- 
fluidity the relation between superfluidity and Bose condensation was permanently 
in the focus of attention, discussing spin superfluidity one cannot avoid to consider 
the concept of magnon Bose condensation, which is vividly debated nowadays. 
Section W0\ addresses this issue. 



Part I: Spin currents in magnetically ordered systems 



2. Mass supercurrents 

Since the idea of spin superfluidity originated from the analogy with the more 
common concept of mass superfluidity let us shortly summarize the latter. The 
essence of the transition to the superfluid or superconducting state is that below 
the critical temperature the complex order parameter ip = \ip\e^^ , which has a 
meaning of the wave function of the bosons or the fermion Cooper pairs, emerges 
as an additional macroscopical variable of the liquid. For the sake of simplicity, we 
restrict ourselves to the case of a neutral superfluid at zero temperature putting 
aside the two- fluid theory for finite temperatures. Then the theory of superfluidity 
tells that the order parameter ^ determines the particle density n = |^p and the 
velocity of the liquid is given by the standard quantum-mechanical expression 

• ^ -X'il:*V^-^V^*) = —V^. (1) 



2m\tp\'^ m 

Thus the velocity is a gradient of a scalar, and any flow is potential. Since the 
phase and the particle number are a pair of canonically conjugate variables, one 
can write down the Hamilton equations for the pair of the canonically conjugated 
variables "phase - density" : 

dip 58 dn 6£ , , 

dt 6n^ dt hSif 

Here £ = J d^RE is the total liquid energy, whereas E is the energy density, and 
5£/5n and 6£/6ip are functional derivatives of the total energy: 

/^, (3) 
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Figure 1. Phase (inplane rotation angle) variation at the presence of mass 
(spin) supercurrents. a) Oscillations in a sound (spin) wave), b) Stationary 
mass (spin) supercurrent. 



In these expressions fi is the chemical potential, 

dE 



nv 



hdVip 



(5) 



is the particle current, and the dependence of the energy on the density gradient 
was ignored. Eventually the Hamilton equations are reduced to the equations of 
hydrodynamics for an ideal liquid: 



m- 



dv 
~dt 



-V//, 



(6) 



dn 
dt 
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A crucial property of the system is the gauge invariance: the energy does not 
depend on the phase directly {dE/dip = 0) but only on its gradient. According to 
Noether's theorem this must lead to the conservation law for a conjugate variable, 
the total number of particles. The conservation law manifests itself in the conti- 
nuity equation ([7]), which contains the particle supercurrent. The prefix "super" 
stresses that this current is not connected with dissipation. It is derived from the 
Hamiltonian or the Lagrangian but not from the dissipation function. In contrast to 
the diffusion current proportional to the density gradient, the supercurrent is pro- 
portional to the phase gradient. Therefore it appears only in a coherent state with 
broken gauge invariance. The equations of superfluid hydrodynamics can be derived 
from the Gross-Pitaevskii equation for a weakly non- ideal Bose-gas. However, they 
are much more general than this model. They can be formulated from the most gen- 
eral principles of symmetry and conservation laws. Indeed, deriving the two-fluid 
theory of superfluidity Landau did not use the concept of Bose-condensation. 

An elementary collective mode of the ideal liquid is a sound wave. In a sound 
wave the phase varies in space, i.e., the wave is accompanied by mass supercur- 
rents (figure flk) . An amplitude of the time and space dependent phase variation 
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is small, and currents transport mass on distances of the order of the wavelength. 
A really superfluid transport on macroscopic distances is related with stationary 
solutions of the hydrodynamic equations corresponding to finite constant currents 
with constant nonzero phase gradients (current states). In the current state the 
phase rotates through a large number of full 27r-rotations along streamlines of the 
current (figure [ija). 

The crucial point of the superfiuidity concept is why the supercurrent is a persis- 
tent current, which does not decay despite it is not the ground state of the system 
and has a larger energy. The first explanation of the supercurrent stability was 
given on the basis of the well known Landau criterion [16j . According to this cri- 
terion, the current state is stable as far as any quasiparticle of the Bose-liquid in 
the laboratory frame has a positive energy and therefore its creation requires an 
energy input. Let us suppose that elementary quasiparticles of the Bose-liquid at 
rest have an energy spectrum e(p) where p is the quasiparticle momentum. If the 
Bose-liquid moves with the velocity v the quasiparticle energy in the laboratory 
frame is e{p) +pv. The energy cannot be negative (which would mean instability) 
if 

e(p) ^_. 

V < vl = mm . (8) 

P 

In superfiuid ''He the Landau critical velocity vl is determined by the roton part 
of the spectrum. But in this review we focus on the long-wavelength collective 
excitations, which are phonons with the spectrum e = UsP- Then according to 
equation (|8| the supercurrent cannot be stable if the velocity v exceeds the sound 
velocity Ug. 

As far as one wants to check the Landau criterion for long-wavelength collective 
modes like sound waves, it is not necessary to solve a dynamical problem looking 
for the spectrum of phonons. It is enough to estimate the energy of possible static 
fluctuations in the stationary current state with particle density no and velocity vq. 
Let us write down the energy of the current state taking into account possible local 
fluctuations of the particle density, n' = n — uq, and of the velocity, v' = v — vq, 
up to the terms of the second order: 



£= d^R 



^o(no)n + -^^^ + 2 

Here //q = dEo/dn and Eq are the chemical potential and the energy density of 
the liquid at rest. One may neglect terms of the first order with respect to the 
density fluctuation n' and the velocity v' since we look for an energy extremum 
at fixed averaged density and velocity and the first-order term must vanish after 
integration. Using the thermodynamic relation mu1/n = dfio/dn = d'^EQ/dn? and 
omitting the subscript in uq and vq, one obtains 



£^ d^R 



1 12 II 
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(9) 



The quadratic form under the integral is positive definite, i.e., the current state 
corresponds to the energy minimum, if the condition u^ > v^ is satisfied. This 
condition is identical to the Landau criterion equation ([s]) for the phonon spectrum 
e = Usp. 

The theory of superfluidity tells that the Landau criterion is a necessary but 
not sufficient condition for current metastability. The Landau criterion checks only 
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Figure 2. Mass and spin vortices, a) Mass vortex or spin vortex in an easy- 
plane ferromagnet without inplane anisotropy. b) Spin vortex at small spin 
currents ((Vi/?) <§; l/l) for four-fold inplane symmetry. The vortex line is a 
confluence of four 90° domain walls (solid lines). 



small deviations from the current state. Meanwhile the current state can be de- 
stroyed via large perturbations of the current state. In superfluids these large per- 
turbations are vortices. In the current state the phase rotates along the current 
direction. The current can relax if one can remove one 27r-turn of the phase. This 
requires that a singular vortex line crossed or "cut" the channel cross-section. The 
process is called "phase slip" . 

If the vortex axis (vortex line) coincides with the z axis, the phase gradient 
around the vortex line is given by 



Vip^ 



[z X r] 



(10) 



where r is the position vector in the xy plane. The phase changes by 27r around 
the vortex line (figure l2k) . Creation of the vortex requires some energy. The vortex 
energy per unit length (line tension) is determined by the kinetic (gradient) energy: 



a r = m — 



2?7l 



m 



(11) 



where the upper cut-off r^ is determined by geometry. For example, for the vortex 
shown in figure [2^ it is the distance of the vortex line from a sample border. The 
lower cut-off re is the vortex-core radius. It determines the distance r at which the 
phase gradient is so high that the hydrodynamic expression for the energy becomes 
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Figure 3. Topology of the uniform mass current and the vortex states, a) The 
current state in a torus maps onto the circumference \ip\ = |-(/>o| = const in the 
complex tp - plane, where ipo is the equilibrium order parameter wave function 
of the uniform state, b) The vortex state maps onto the circle \ip\ < |i/'o|- 



invalid. A good estimation for re is Tc ~ n/us, where k = h/m is the circulation 
quantum of the velocity. Inside the core the modulus of the order parameter goes 
down to zero eliminating the singularity in the kinetic energy at the vortex axis. 
For the weakly non-ideal Bose-gas this estimation yields the coherence length. 

Now suppose that a vortex appears in the current state with the constant gradient 
V(/?o: The phase gradients induced by the vortex are superimposed on the constant 
phase gradient related to the current: V93 = V990 + V(/?^. The total gradient 



energy includes that of the current, the vortex energy given by equation (11), 
and the energy from the cross terms of the two gradient fields. Only the last two 
contributions are connected with the vortex, and their sum determines the energy 
of the vortex in the current state: 

g="^Lln^-'^5Vv.o, (12) 

m re m 

where L is the length of the vortex line and S is the area of the cut, at which 
the phase jumps by 27r. For the 2D case shown in figure ^ (a straight vortex in 
a slab of thickness L normal to the picture plane) S = Lrm- One can see that 
vortex motion across the channel (growth of r^) is impeded by the barrier, which 
is determined by variation of the energy e with respect to r^. The peak of the 
barrier corresponds to r^ = l/2V(/3o- The height of the barrier is 

irh'^n 1 

em ~ im— - — . (13) 

m rcVifo 

Thus the barrier disappears at gradients V(/9o ~ l/'''c, which are of the same order 
as the critical gradient determined from the Landau criterion. In the 3D geometry 
the phase slip is realized with expansion of vortex rings. For the ring of radius R 
the vortex-length and the area of the cut are L = 2'kR and 5 = vri?^ respectively, 
and the barrier disappears at the same critical gradient ~ 1/rc as in the 2D case. 
The barriers stabilizing metastable current states are connected with topology 
of the order parameter space. In a superfluid the order parameter is a complex 
wave function ^{r). At the equilibrium tp = |'0o|e*''') where the modulus iV'ol is a 
constant determined by minimization of the energy and the phase c^ is a degen- 
eration parameter since the energy does not depend on 99. Any current state in 
a closed annular channel (torus) with the phase change 27rn around the channel 
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maps onto a circumference \tp\ = iV'ol in the complex plane (figure^) winding the 
circumference n times. It is evident that it is impossible to change n keeping the 
path on the circumference \ip\ = iV'ol ^ill the time. Thus n is a topological charge. 
One can change it (removing, e.g., one winding around the circumference) only 
by leaving the circumference (the equilibrium order parameter space in the case). 
This should cost energy, which is spent on creation of a vortex. Figure [Sja shows 
mapping of the vortex state onto the circle \ip\ < IV'ol- 

Without such topological barriers superfluidity is ruled out. However, barriers 
do not automatically provide the life-time of currents long enough. In practice, 
dissipation via phase slips is possible even in the presence of barriers due to thermal 
fluctuations or quantum tunneling. Here and later on we address only "ideal" 
critical currents (the upper bound for critical currents) at which barriers disappear 
leaving "practical" critical currents beyond the scope of the present review. 



3. Phenomenology of magnetically ordered systems and spin currents 

The main interaction responsible for magnetic order is exchange interaction, which 
is invariant with respect to rotations of the whole spin system. Then according 
to Noether's theorem the total spin must be conserved. For ferromagnets where 
the order parameter is the spontaneous magnetization M, this means that the 
exchange energy can depend on the absolute value of M but not on its direction. 
Other contributions to the free energy (anisotropy energy or dipole-dipole inter- 
action) are related with spin-orbit interaction, which does not conserve the total 
spin. But these interactions are relativistically small, i.e., governed by the small 
relativistic parameter v/c, where t; is a typical electron velocity and c is the speed 
of light. The spin-orbit interaction does depend on M direction, but because of its 
weakness cannot affect the absolute value M in slow dynamics. This is a crucial 
point in the phenomenological theory of magnetism of Landau and Lifshitz [17|. 
which determines the form of the equation of motion for ferromagnet magnetization 
known as the Landau-Lifshitz equation [18j : 

dM 

^ = j[H,jfXM], (14) 

where 7 is the gyromagnetic ratio between the magnetic and mechanical moment 
(M = —jS). The effective magnetic field is determined by the functional derivative 
of the total free energy J~ = f d^RF with density F: 

H.„ = -^. (15) 

According to the Landau-Lifshitz equation, the absolute value M of the magnetiza- 
tion cannot vary. The evolution of M is a precession around the effective magnetic 
field Heff . 

At first let us discuss exchange approximation, in which relativistic effects are 
ignored and the conservation of total spin is not violated. In this approximation 
the free energy density is 

The first exchange-energy term Fq{M), being the largest term, is crucial for deter- 
mination of the equilibrium value of M. But after determination of M it can be 
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ignored as an inessential large constant. Indeed, its contribution to the effective field 



in Landau-Lifshitz equation (14) does not produce any effect: the contribution is 
parallel to M and vanishes in the vector product. In the absence of external fields, 
which break invariance with respect to rotations in the spin space, the Landau- 
Lifshitz equation reduces to the continuity equations for components of the spin 
density S = — iW/7: 

where 

dF 






M X 



dVjM 



a[M X VjM]i = -aSikiMkVjMi (18) 



is the jth component of the spin current transporting the ith component of spin. 
Thus in an isotropic ferromagnet all three components of spin are conserved. 

The Landau-Lifshitz equation has plane-wave solutions describing spatially 
nonuniform precession of the magnetization M = Mq + m around the ground- 
state magnetization Mq: m oc e**"""*"^*. Here the magnetization deviation m is 
small and normal to Mq. Linearizing with respect to m, one obtains spin waves 
with the spectrum 

u = jaMok'^. (19) 

In an isotropic ferromagnet spin waves at A: 7^ are accompanied by spin currents, 
but superfluid spin transport is impossible as will be clear from section |4} 

Next we shall consider the case when spin-rotational invariance is partially bro- 
ken, and there is uniaxial crystal magnetic anisotropy given by the third term in 
the phenomenological free energy: 

a M^ 

F = Fo{M) + -ViM-ViM + EA^. (20) 

If the anisotropy energy E^ is positive, it is the "easy plane" anisotropy, which 
keeps the spontaneous magnetization Mq in the xy plane (the continuous limit 
of the XY model). In this model the z component of spin is conserved, because 
invariance with respect to rotations in the easy plane remains unbroken. Since the 
absolute value of magnetization is fixed, the vector M of the magnetization is 
fully determined by the angle (/? showing the direction of M in the easy plane xy 
{Mx = Mcosc/7, My = M simp) and by the z component of the magnetization rriz- 
We use the notation niz instead of M^ in order to emphasize that ruz is a small 
dynamic correction to the magnetization, which is absent at the equilibrium. In 
the new variables the free energy is 



-F= / d^RF = / d^R 



2x 2 



(21) 



The constant A = aM"^ is stiffness of the spin system determined by exchange 
interaction, and the magnetic susceptibility x = M'^ /Ea along the z axis is de- 
termined by the uniaxial anisotropy energy Ea keeping the magnetization in the 
plane. The Landau-Lifshitz equation reduces to the Hamilton equations for a pair of 
canonically conjugate continuous variables "angle-angular momentum" (analogous 
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to the canonically conjugate pair "coordinate-niomentum" ) : 



^ = -7— = -7— (22) 

dt Sruz druz ' 



Idruz 6T OF „ OF 

-^ = ^ = 7^-^-7^^^' 23 

7 at dip dip ovip 

where functional derivatives on the right-hand sides are taken from the free energy 



F given by equation (21). Using the expressions for functional derivatives one can 



write the Hamilton equations as 



^._7- (24) 



''^"^^+V.J^ = 0, (25) 



7 dt 



where 



dF 

'dVip 



-AVif (26) 



is the spin current. 

There is an evident analogy of equations (24) and (25) with the hydrodynamic 
equations ^ and ([T]) for an ideal liquid, equation ( 25 ) being the continuity equa- 
tion for spin. This analogy was exploited by Halperin and Hohenberg [T9] in their 
hydrodynamic theory of spin waves. In contrast to the isotropic ferromagnet with 
the quadratic spin- wave spectrum, the spin wave in the easy-plane ferromagnet has 
a sound-like spectrum as in a superfluid: uj = Csk, where the spin-wave velocity 
is Cs = 7yM7x- Halperin and Hohenberg introduced the concept of spin current, 
which appears in a propagating spin wave like a mass supercurrent appears in a 
sound wave (figure Ilk) . This current transports the z component of spin on dis- 
tances of the order of the wavelength. But as well as the mass supercurrent in 
a sound wave, this small oscillating spin current does not lead to superfluid spin 
transport, which this review addresses. Spin superfluid transport on long distances 
is realized in current states with magnetization rotating in the plane through a 
large number of full 27r-rotations as shown in figure [T]3. 

Let us consider now the case of antiferromagnetic order. The simplest model of 
an antiferromagnet is two sublattices with magnetizations Mi and M2- In the 
absence of weak ferromagnetism and external magnetic fields two magnetizations 
Ml = —M2 completely compensate each other without producing any total mag- 
netization m = Ml + M2 ■ However, a small magnetization m does appear due to 
external magnetic fields or dynamical effects. The amplitudes of Mi and M2 and 
their mutual orientation are mostly determined by strong exchange interaction, but 
the latter does not fix the direction of the staggered magnetization L = Mi — M2, 
which is the order parameter of a two-sublattice antiferromagnet. The equations 
of motion for two vectors L and m can be derived from the two Landau-Lifshitz 
equations for Mi and M2 taking into account the exchange interaction between 
two sublattices. But it would be useful to present a more general version of the 
macroscopic phenomenological theory, which is able to describe an antiferromag- 
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netic structure of any complexity [20^ \2T\ . The theory of spin dynamics in su- 
perfluid phases of '^He developed by Leggett and Takagi [22] also belongs to this 
class. Following the same principle "exchange is the strongest interaction" as in 
the Landau-Lifshitz theory, macroscopic theories of this type deal with phenomena 
at scales essentially exceeding microscopic scales (the coherence length in the case 
of ^He), at which the exchange energy establishes the tensor structure of the order 
parameter. This permits to assume that the entire dynamic evolution of the order 
parameter reduces to rotations in the 3D spin state, which cannot change the ex- 
change energy. Then the dynamics of the system is described by three independent 
pairs of canonically conjugated variables "angle-moment" ipi-rrii/^ (i = 1,2,3): 

d(pi 6F 



dt Srrii 

1 drrii 5F 



7- 
trrii 

(27) 



7 dt 6(pi 

Here ipi are the angles of spin rotations around three Cartesian axes (i = x,y,z). 
Apart from spatial dependence of the variables, these equations are similar to the 
equations of motion of a 3D rigid top. In our case the top is an antiferromagnetic 
spin order parameter rigidly fixed by exchange interaction. As in the case of a 
two-sublattice antiferromagnet, magnetization m results from deformation of the 
equilibrium spin structure. The approach is valid as far as this deformation is weak, 
i.e., m is smaller than the characteristic moments of the antiferromagnetic structure 
(staggered magnetization L in the case of a two-sublattice antiferromagnet). Since 
rotation around the vector L has no effect on the state of the system the latter has 
only two degrees of freedom corresponding to two pairs "angle-moment". Then 
the equations become the equations of motion of a rotator. In contrast to the 



spontaneous magnetization M in the Landau-Lifshitz equation (14), the small 
absolute value of the magnetization m is not kept constant. 

Because the group of 3D rotations is non-commutative, the state of the system 
depends on the order, in which rotations around different axes are performed. 
In practice they frequently use the Euler angles (they are introduced in section 
nh. For the most content of Part I (except for section wn, one can choose one 
degree of freedom connected with the conjugate pair c^^-m^, and the problem of 
non-commutativity is absent (further we shall omit the subscript z of the angle 
ifz)- If the energy of the ground state does not depend (or depends weakly as 
discussed in section pi) on the angle ^p, the equations of motions for 93 and niz are 



the same Hamilton equations (22) and(23), which were formulated for an easy 



plane ferromagnet. In the case of a two-sublattice antiferromagnet the angle (/? is 
the angle of the staggered magnetization L in the easy plane. 

The discussion of this section has not made any reference to a concrete micro- 
scopic model of magnetism. Indeed, the approach is general enough and is valid 
for models of magnetism based on the concepts of either localized or itinerant 
electrons. In particular, ferromagnetism of localized electrons is described by the 
Heisenberg model with the Hamiltonian: 

H = -JJ2^^■^i+l^ (28) 

where J > 0, Sj are spins at the sites i, and the summation over j includes only the 
nearest neighbors to the site i. In the continuum limit, when the spin rotates very 



slowly at scales of the intersite distance a, the Hamiltonian ( 28 ) reduces to the free 
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energy (16) in the Landau-Lifshitz theory with the magnetization M = —^{si)/a'^ 
and the stiffness constant a = Ja^/7^. 

The debates on rehabihty of the general phenomenological approach to mag- 
netism are as old as the approach itself. Nearly sixty years ago Herring and Kit- 
tel [23] argued with their opponents that their phenomenological theory of spin 
waves "is not contingent upon the choice of any particular approximate model 
for the ferromagnetic electrons". Interestingly these discussions are still continu- 
ing in connection with the concept of the spin current, which originates from the 
general phenomenological approach. Originally they connected spin supercurrents 
with counterflows of particles with opposite spins, for example, of He atoms in 
the A-phase of ^He [HIS]. Bunkov ^3\ insisted that only a counterflow of particles 
with opposite spins would lead to superfluid spin transport, thus ruling out spin 
superfluidity in materials with magnetic order resulting from exchange interaction 
between localized spins (see p. 93 in his review). However, the spin current does not 
require itinerant electrons for its existence [2] . The presumption that spin transport 
in insulators is impossible is still alive nowadays. According to Shi et al. L24J, it is 
a critical flaw of spin-current definition if it predicts spin currents in insulators. 



4. Stability of spin-current states 

For the sake of simplicity further we focus on current states in an easy-plane fer- 
romagnet, though the analysis can be easily generalized to other magnetically or- 
dered systems discussed in the previous section. In the current state the sponta- 
neous magnetization M(r) rotates in the easy plane through a large number of 
full 27r-rotations when the position vector r is varying along the direction of spin 
current (figure lip). The spin-current state is metastable if it corresponds to a local 
minimum of the free energy, i.e., any transition to nearby states would require an 
increase of energy. This condition is an analog of the Landau criterion for mass 
supercurrents discussed in section [2] In order to check current metastability, one 
should estimate the energy of possible small static fluctuations around the sta- 
tionary current state. For this estimation, one should take into account that the 
stiffness constant A is proportional to the squared inplane component of the spon- 
taneous magnetization Mj_ = Mq — ml , and in the presence of large angle gradients 
A must be replaced with A[l — tti^/Mq). So the free energy is 



T = / d^R 



d^R 



ml A{l-ml/Mi){V^r 
2x 2 

mlEA-A{V^f A{V^f^ 



Mi 



+ 



(29) 



One can see that if V93 exceeds ^J M'^/xA = ^J EaI A the current state is unstable 
with respect to the exit of iWo from the easy plane. This is the Landau criterion 
for the stability of the spin current. 

Like in superfluids, stability of current states is connected with topology of the 
order parameter space. For ferromagnets the order parameter is the magnetization 
vector NL. For isotropic ferromagnets the space of degenerated equilibrium states 
is a sphere \M\ = |-/Vfo|, whereas for an easy-plane ferromagnet this space reduced 
to an equatorial circumference on this sphere (figure [4k). Thus the order parameter 
space for an easy-plane ferromagnet is topologically equivalent to that space for su- 
perfluids (the circumference on the complex plane shown in figure p|. Spin-current 
states are stable because they belong to the topological classes different from the 
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Figure 4. Topology of the uniform spin current and the spin vortex states, 
a) The current state in a torus maps onto the equatorial circumference of 
the order parameter sphere \M\ = const, b) Isotropic fcrromagnet: continu- 
ous deformation reduces a circumference (current state) to a point (uniform 
current-free state), c) The spin vortex state maps onto either an upper or a 
lower half of the sphere \M\ = const. 



class of the uniform ground state and cannot be reduced to the latter by continuous 
deformation of the path. In contrast, for an isotropic ferromagnet the path around 
the equatorial circumference can be continuously transformed to a point on the 
sphere as shown in figure [Ija. In this process the energy monotonously decreases, 
and topological barriers are absent. Topology of an easy-axis (anti)ferromagnet 
also does not allow stable spin-current states. 

The connection of superfluidity-like phenomena with topology of the order pa- 
rameter space is universal and not restricted with the examples of mass and spin 
superfluidity considered here. The same arguments support possibility of exciton 
superfluidity, which was discussed even earlier than spin superfluidity (see the in- 
troductory section Qj). Though the whole problem of superfluid exciton transport 
is far from its resolution, in some special case experimental evidences of this trans- 
port has already been reported. Kellogg et al. [25] observed vanishing resistance 
in double quantum Hall layers, which was interpreted as a consequence of Bose 
condensation of interlayer excitons (or pseudospin ferromagnetism) . 

As well as in the theory of mass superfluidity, the Landau criterion is a nec- 
essary but not sufficient condition for current metastability. One should also to 
check stability with respect to large perturbations, which are magnetic vortices. 
The magnetic vortices were well known in magnetism. Bloch lines in ferromagnetic 
domain walls are an example of them |26j . In the spin-current state the magneti- 
zation M traces a spiral at moving along the current direction. The spin current 
can relax if one can remove one turn of the spiral. This requires that a singular 
line (magnetic vortex) crossed or "cut" the channel cross-section [2i i3j as shown in 
figure [2|3. 

The structure of the magnetic vortex outside the vortex core is the same as of 



March 3, 2010 17:4 Advances in Physics SpinRev 



14 E. B. Sonin 



the mass superfluid vortex given by equation (10). Correspondingly, the magnetic 



vortex energy is determined by the expression similar to equation (11): 



dV^i^^^=vr^ln^, (30) 

2 re 

where the upper cut-off Vm depends on geometry. However, the radius Vc and the 
structure of the magnetic vortex core are determined differently from the mass 
vortex. In a magnetic system the order parameter must not vanish at the vortex 
axis since there is a more effective way to eliminate the singularity in the gradient 
energy: an excursion of the spontaneous magnetization out of the easy plane xy. 
This would require an increase of the uniaxial anisotropy energy, which keeps M 
in the plane, but normally this energy is much less than the exchange energy, 
which keeps the order-parameter amplitude M constant. Finally the core size re 
is determined as a distance at which the uniaxial anisotropy energy density Ea is 
in balance with the gradient energy A{'Vip)'^ ~ ^/''"c- This yields Vc ~ \J AJEa- 
Figure |4]3 shows mapping of the spin vortex state onto the order parameter space. 
In contrast to superfluid vortices mapping onto a plane circle, the spin vortex state 
can map onto one of two halves of the sphere \M\ = const. Thus a magnetic (spin) 
vortex has an additional topological charge having two values ±1 [27]. 

The energy of the spin-current state with a vortex and the energy of the barrier, 
which blocks the phase slip, i.e., the decay of the current, are determined similarly 



to the case of mass superfluidity [see equations (12) and (13)]: 



e = TTLAln—-2TTASV^o, (31) 



e^^vrLAln— ^, (32) 

where L is the length of the vortex line and S is the area of the cut, at which the 
angle jumps by 27r. Thus the barrier disappears at gradients Vipo ~ l/^'o which are 
of the same order as the critical gradient determined from the Landau criterion. 
This is a typical situation in the superfluidity theory. But sometimes the situation 



is more complicated as we shall see in section 7.2 



5. Spin currents without spin conservation law 

Though processes violating the conservation law for the total spin are relativisti- 
cally weak, their effect is of principal importance and in no case can be ignored. The 
attention to superfluid transport in the absence of conservation law was attracted 
first in connection with discussions of superfluidity of electron-hole pairs. The num- 
ber of electron-hole pairs can vary due to interband transitions. As was shown by 
Guseinov and Keldysh |28j . interband transitions lift the degeneracy with respect 
to the phase of the "pair Bose-condensate" and make the existence of spatially 
homogeneous stationary current states impossible. On the basis of it Guseinov and 
Keldysh concluded that there is no analogy with superfluidity. This phenomenon 
was called "fixation of phase" . However some time later it was demonstrated |29j 
that phase fixation does not rule out existence inhoniogeneous stationary current 
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states, which admit some analogy with superfluid current state^ This analysis 
was extended on spin currents [21 13] . 

In the spin system the role of the phase is played by the angle of the magnetization 
M in the easy plane, and the degeneracy with respect to the angle is lifted by 
magnetic anisotropy in the plane. Adding the n-fold inplane anisotropy energy to 



the total free energy (21) the latter can be written as 
T = I (TH { -^ + 



m% Aiy^pY 



2x 



+ K\V — cos(nc/7)] 



(33) 



Then the spin continuity equation (25) becomes 
1 dmz 



7 dt 



V • J^ + nK'smirup) = -A 



'y2 _ sin(rav?) 



where 



I' 



A 



(34) 



(35) 



Excluding m^ from equations (24) and (34) one obtains the sine Gordon equation 
for the angle (p: 



av 
at2 



■^2 _ s\n{nip) 



0, 



(36) 



where Cg = yJ^A/x is the spin-wave velocity. According to this equation, the 
inplane anisotropy leads to a gap in the spin-wave spectrum: 



UJ 



-"2 h Cg/e . 



(37) 



There are one-dimensional solutions 93(2; — vt) of the sine Gordon equation with 
non-zero average (Vc/9), which correspond to a periodic lattice of solitons (domain 
walls) of the width ~ I = 1^1 — v"^ /c^ with the period xq = 27r/n{V(p) moving 
with the velocity v. The function inverse to ^{x — vt) is 



X — vt 



"[ 



dip' 



2 J ^K — cos TUp' 

where the constant k > 1 is determined by the equation 



27r 



xo 



^[ 



27r/n 



d^' 



niyp) V 2 Jo a/^ — cos nip' 
The free energy of the soliton lattice is given by 



F = A 



1 - K {Vp) [^ /"^-Z" ^ ^ 

H ^-v/— / \Jk — cosmpap 



nP 



■nl 



(38) 



(39) 



(40) 



^Similar conclusions have been done with respect to possibiUty of supercurrents in systems with spatially 
separated electrons and holes 1301 131| . 
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Figure 5. The nonuniform spin-current states with (Vi/s) -C l/l and (V</3) 2> 
l/l. 



It is possible to develop the hydrodynamic theory of the soliton lattice in the 
terms of local density and velocity of solitons [32], which is able to describe defor- 
mations of the lattice slow in space and time. Here we focus on stationary current 
states when drriz/dt = (f = 0). At small average twisting of the spontaneous 
magnetization (Vip) <C l/l the structure constitutes domains that correspond to 
the n equivalent easiest directions in the easy plane. In this limit (k — )• 1) the free 
energy density is the product of the energy of an isolated domain wall and the 
density of domain walls n(V(^)/27r: 



F = A 



4(Vy) 
TT\/nl 



(41) 



Spin currents (gradients) inside domains are negligible but there are essential spin 
currents inside domain walls where Vip ~ l/l. This hardly reminds genuine su- 
perfluid transport on macroscopical scales: spin is transported over distances on 
the order of the domain- wall width /. With increasing (Vt/?) the density of domain 
walls grows, and at {V(p) ^ l/l they coalesce while for a displacement along the 
direction of the gradient (V93), the end point of the vector M describes, a line close 
to a helix. The nonuniform states with (Vip) <C l/l and (Vip) S> l/l are shown in 
figure [5} Thus the processes violating spin conservation law are not important for 
large deformations (gradients) of the spin structure. This means that the analogy 
of these deformed states with current states in superfluids makes sense. 

Studying stability of nonuniform current states it is possible to ignore the inplane 
anisotropy only for large spin currents when V99 ^ l/l. Let us consider the opposite 
limit of (V99) <^ l/l when the spin structure reduces to a chain of domain walls. The 
relaxation of the spin current, which is proportional to the wall density, requires 
that some domain walls vanish from the channel. This process is illustrated in figure 
pja for the four-fold inplane symmetry (n = 4). When a magnetic vortex appears, n 
domain walls finish not at the wall but at the vortex line, around which the angle 
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ip changes by 2tt. The 2tt angle jump occurs at the cut restricted by the vortex 
hne. The n domain wahs disappear via motion of the vortex line across the channel 
cross section. In the course of this process, the change of the energy consists of the 
vortex-line energy, which is proportional to the line length L, and of a decrease 
of the surface energy of the n walls themselves proportional to the cut area S. 
The latter contribution is determined by the product of the free energy density 



(41) and the volume 27rS'/(Vc/9). Taking these two contributions into account, the 



energy during the process of annihilation of n walls is 



i = nLAln'^-~. (42) 



Comparing it with the energy given by equation (31) one sees that the gradient 



V(po is replaced by the maximum gradient ~ l/l inside the domain wall. Corre- 



spondingly for the 2D case shown in figure pb the expression ( 32 ) for the barrier 
energy must be replaced by 

e^^^L^ln-^^L^ln^, (43) 

r,, 2 K 



where the two lengths re ~ ^ AJEa and / ~ ^ AjK are determined by the uniaxial 
and the inplane anisotropy energies Ep^ and K. Thus large barriers stabilizing spin- 
current states are possible only if the condition E^ ^ i^ is satisfied. This conclusion 
[21 O [11] was recently confirmed by the analysis of Konig et al. [33]. 

An important difference with conventional mass superfiuidity is that in conven- 
tional superfiuidity the barrier, which suppresses supercurrent relaxation, grows 
unrestrictedly when the gradient V93 decreases. In contrast, in spin superfluidity 
the barrier growth stops when the gradient reaches the values of the order 1// (in- 
verse width of the domain wall). Since the current relaxation time exponentially 
depends on the barrier (whether the barrier is overcome due to thermal fluctuations 
or via quantum tunneling) the life time of the current state in conventional super- 
fluidity diverges when the velocity (phase gradient) decreases. In contrast, the life 
time of the spin current can be exponentially large but always finite. This provides 
an ammunition for rigorists, who are not ready to accept the concept "spin super- 
fluidity" (or superfluidity of any non-conserved quantity) in principle. In principle, 
one could agree with them. But in practice, whatever we call it, "non-ideal su- 
perfluidity" or "quasi-superfluidity" , some consequences should outcome from the 
fact of the existence of topological barriers suppressing relaxation of spin-current 
states. A key point is whether these consequences are observable. This is the topic 
of the next section. 



6. Is superfluid spin transport "real"? 

From early days of discussions on spin supercurrents and up to now there are ar- 
guments on whether the spin supercurrent can result in "real" transport of spin. 
Partially this is a semantic problem: One must carefully define what "real" trans- 
port really means. Let us suppose that one has a usual superfiuid mass persistent 
current in a ring geometry. Nobody doubts that real mass transport occurs in this 
case, but how can one notice it in the experiment? In any part of the ring channel 
there is no accumulation (increase or decrease) of the mass. Of course, one can 
detect gyroscopic effects related with persistent currents, but it is an indirect evi- 
dence. What may be a direct evidence? One could suggest a Gedanken Experiment 
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Supetflow of angular momentum 



Figure 6. Mechanical analogue of a persistent current: A twisted elastic rod 
bent into a closed ring. There is a persistent angular-momentum flux around 
the ring. 



in which the ring channel is suddenly closed in some place. In the wake of it one 
can observe that the mass increases on one side from the closure and decreases on 
the other side. This would be a real transport if one required a demonstration of 
mass accumulation as a proof of it. Accepting this definition of transport reality 
one can notice real transport only in a non-equilibrium process, when the trans- 
ported quantity decreases in some place and increases in another . Naturally one 
can discard these semantic exercises as irrelevant for practice, but only as far as 
they refer to mass currents. In the case of spin currents in the past and nowadays 
spin accumulation sometimes is considered as a necessary proof of real spin trans- 
port. Therefore, in old publications on spin superfluidity [21 [3l [11] much attention 
was paid to possible experimental demonstration of spin transport from one place 
to another. 

Before starting discussion of possible spin-transport demonstration it is useful to 
consider a mechanical analogue of superfluid mass or spin supercurrent III] . Let us 
twist a long elastic rod so that a twisting angle at one end of the rod with respect 
to an opposite end reaches values many times 2tt. Bending the rod into a ring 
and connecting the ends rigidly, one obtains a ring with a circulating persistent 
angular- momentum flux (figure ^. The intensity of the flux is proportional to 
the gradient of twisting angle, which plays the role of the phase gradient in the 
mass supercurrent or the spin-rotation-angle gradient in the spin supercurrent. 
The analogy with spin current is especially close because spin is also a part of the 
angular momentum. The deformed state of the ring is not the ground state of the 
ring, but it cannot relax to the ground state via any elastic process, because it 
is topologically stable. The only way to relieve the strain inside the rod is plastic 
displacements. This means that dislocations must move across rod cross-sections. 
The role of dislocations in the twisted rod is the same as the role of vortices in 
the mass or spin current states: In both of the cases some critical deformation 
(gradient) is required to switch the process on. There are various ways to detect 
deformations or strains in an elastically deformed body. Similarly, it is certainly 
possible, at least in principle, to notice deformation (angle gradient) of the spin 
structure in the spin-current state. It would be a legitimate evidence of the spin 
current, not less legitimate than a magnetic field measured around the ring as an 
evidence of the persistent charge current in the ring. 

Of course, it is not obligatory to discuss the twisted rod in terms of angular- 
momentum flux. One can describe it only in terms of deformations, stresses, and 
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Figure 7. Spin injection to a spin-nonsupcrfluid and a spin-supcrfluid 
medium. 



elastic stiffness. So we must have in mind that there are two languages, or descrip- 
tions of the same physical phenomenon. A choice of one of them is a matter of taste 
and tradition. For example, in order to describe the transfer of momentum they 
use the momentum-flux tensor ( "flux" , or "current" language) in hydrodynamics, 
while in the elasticity theory they prefer to call the same tensor as stress tensor. In 
principle one can avoid the term "superfluidity" and speak only about the "phase 
stiffness" even in the case of mass supercurrents. 

Let us return to possible demonstration of "real" spin transport. Suppose that 
spin is injected into a sample at the sample boundary x = (figureul). The injection 
can be realized practically either with an injection of a spin-polarized current (for 
the sake of simplicity we put aside the problem what happens with charge in this 
case), or with pumping the spin with a circularly polarized microwave irradiation. 
If the medium at x > cannot support superfluid spin transport, the only way of 
spin propagation is spin diffusion described by the equations 



dniz 

Ik 



7V-J^ + 






0, J' 



D.Vm, 



(44) 



where Ds is the spin-diffusion coefficient and Ti is the time characterizing the Bloch 
longitudinal relaxation, which violates the spin-conservation law. In the stationary 
case drriz/dt = 0, and both the spin current and the nonequilibrium magnetization 
rriz exponentially decay inside the sample: JJ oc ttt-^ fx e~^' " , where Lg = y/DgTi 
is the spin-diffusion length. So no spin can reach the other boundary x = L oi the 
sample provided L ^ Lg. 

Now let us suppose that the medium at < x < L is magnetically ordered and 
can support superfluid spin transport. If the injection is so weak that the angle gra- 
dient Vy? is much less than l/l, the perturbation of the medium by the injection can 
penetrate at the length not longer than the domain- wall width /. So the injection 
pushes a piece of a domain wall into the sample. The spin current exponentially 
decays like in the medium without superfluidity. If the injection is so strong that 
the angle gradient Vc^ exceeds its maximum value Yj ^/nl in a center of an isolated 
domain wall, continuity of the spin current on the boundary requires appearance 
of a soliton lattice with a period of the order or smaller than I. This means that 
the spin current can reach the other boundary x = L. In the thermodynamic limit 
L — )• oo one may expect a stationary soliton lattice with ruz = 0. However, at the 
boundary x = L the spin current must be injected into the medium without spin 
superfluidity. This is impossible without a finite ttt,^, and the finite ruz means that 
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the soliton lattice is moving. In the hmit V99 ^ l/l the sohton lattice is rather 
dense and one may neglect periodical spin-current modulation caused by inplane 
anisotropy. Then spin transport is described by the equations 

-^ = -7—, (45) 

at X 



"^-tV.J' + ^-O. (46) 

with the boundary conditions for the supercurrent J^{0) = Jq at x = and 
J^{L) = —fmz{L) at X = L. The current Jq in the first condition is the spin- 
injection current, while the second boundary condition takes into account that 
the medium at x > L is not spin-superfluid and spin injection there is possible 
only if some non-equilibrium magnetization mz{L) is present. The coefficient / can 
be found by solving the spin-diffusion equations in the medium at x > L [3]. It 
also depends on properties of the contact at x = L. While the inplane anisotropy 
violating the spin conservation (phase fixation) was neglected, one cannot neglect 
irreversible dissipative processes, which also violate the spin-conservation law. The 
simplest example of such a process is the longitudinal spin relaxation characterized 
by time Ti. 



The stationary solution of equations (45) and (46) is 



Though the solution is stationary in the sense that druz/dt = 0, but dip/dt / 0. 
We consider a non-equilibrium process (otherwise spin accumulation is impossible) , 
which is accompanied by the precession of M in the easy plane. But the process 
is stationary only if the precession angular velocity is constant in space. The con- 
dition niz =const, which results from it, is similar to the condition of constant 
chemical potential in superfiuids or electrochemical potential in superconductors 
in stationary processes. If this condition were not satisfied, there would be steady 



growth of the angle twist as is evident from equation ( 45 ) . As already mentioned 
above, the nonzero rUz means that the soliton lattice is moving. In our case the 
soliton velocity is rather slow since it is inversely proportional to L\ v = c^Ti/L. 

One can see that irreversible loss of spin is a more serious obstacle for super- 
fluid spin transport than coherent phase fixation, to which most of attention was 
attracted in the literature. Because of spin relaxation, the spin current inevitably 
decreases while moving away from the injection point, in contrast to constant su- 
perfluid mass currents. However, in a spin-superfluid medium this decrease is linear 
and therefore less destructive than exponential decay of currents in non-superfluid 
media. So one have a good chance to notice spin accumulation in the medium at 
X > L rather distant from the place of original spin injection. This justifies using 
the term "superfluid" . 

In the presented analysis we assumed that at the boundary the entire spin- 
injection current is immediately transformed into a supercurrent. Actually spin 
injection can also generate the diffusion current close to the boundary. However, at 
some distance from the boundary the diffusion current inevitably transforms into 
a supercurrent ^ . If this distance (healing length) is much shorter than the size L 
of the sample our boundary condition at x = is fully justified. Similar effects take 
place at contacts "normal metal - superconductor" : The current from the normal 
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metal to the superconductor is completely transformed into the supercurrent at 
some finite distance from the contact. 

Spin injection is not the only method of generation of spin currents. One can 
generate spin currents by a rotating inplane magnetic field, which is applied to 
one end of the sample and is strong enough to orient the magnetization parallel 
to it. Because of the stiffness of the spin system, the spin rotation at one end is 
transmitted to the other end of the sample, which is not subject to the direct effect 
of the rotating magnetic field. Transmission of the torque through the sample is 
spin current. Since the rotating field is acting on the phase (angle) of the order 
parameter, it can be called coherent method, in contrast to the incoherent method 
of spin injection. The coherent method of spin-current generation has no analog in 
superfluids and superconductors, since in the latter cases there is no field linked 
to the phase of the order parameter. Referring to the set up shown in figure [7] 
with spin injection replaced by rotating magnetic field, in the coherent method the 
magnetization ruz is fixed by the frequency of the rotating field. In this case there 
is no threshold for spin-current generation, and the spin current appears whatever 
small the frequency could be. But if the frequency (and niz proportional to it) is low, 
the spin transmission is realized via generation of a chain of well separated solitons 
(domain walls), which propagate to the other end of the sample. Thus, a "moving 
soliton lattice" is another synonym for spin superfluid transport. Long-distance 
propagation of solitons through a slab of the A phase of superfluid ^He generated 
by a pulse of a radio-frequency magnetic field has already experimentally realized 
by Bartolac et al. ^34j- This experiment was discussed in terms of spin transport 
in reference 1321. 



7. Spin-precession superfluidity in superfluid ^He-B 

7.1. Stationary uniform precession in ^He-B 

Now we focus on the experimental and theoretical investigations of superfluid spin 
transport in the B phase of superfluid ^He. The spin superfluidity in the B phase 
has several important features, which distinguish it from the spin superfluidity dis- 
cussed previously in this review. First, in contrast to what was considered earlier, 
observed spin-current states in the B phase are dynamical nonlinear states very 
far from the equilibrium, which require for their support permanent pumping of 
energy. Thus dissipation is always present, and speaking about "superfluidity", i.e., 
"dissipationless" spin transport, we have in mind the absence of additional dissipa- 
tion connected with the spin current itself. Second, while the previous discussion 
dealt with the transport of a single spin component (z-component), in the B phase 
spin vector performs a more complicated 3D rotation and the spin current refers 
to the transport of some combination of spin components. This combination may 
be called "precession moment" because it is a canonical conjugate of the preces- 
sion rotation angle (precession phase) rather the rotation angle of genuine spin 
in the spin space. So one should discern two types of spin superfluid transport: 
transport of spin precession and transport of spin |35) . In the experiment [H] they 
used slightly nonuniform magnetic fields, and precession took place only inside 
the homogeneously precessing domain (HPD). But for discussion of superfluid spin 
transport it is not so important, and in the following we consider only processes 
inside the HPD ignoring gradients of the magnetic field. 

The spin dynamics of superfluid phases of ^He is described by the theory of 
Leggett and Takagi [22], which is an example of the general phenomenological 
theory of magnetically ordered systems in terms of conjugate canonical variables 
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Figure 8. Euler angles for spin precession in ^He-B. The rigid top presents 
the rigid spin order parameter structure controlled by large exchange energy. 
The rotation of the coordinate frame transforming the axis z to the axis ^ is 
performed around the axis A'^ (line of nodes). 



"angle-moment" , which was shortly discussed in section |3j As well as in studies 
of rotating solid tops, sometimes it is more convenient to describe spin rotations 
via the Euler angles a, /3, and F (figure^. In ^He superfluid spin dynamics these 
angles were used by Fomin [Hllll], who, however, replaced the angle F by the angle 
<1> = a + F. The angle j3 is the precession tipping angle, and a is the precession 
phase determining the direction of the line of nodes A^. The angle <I> characterizes 
the resultant rotation of the order parameter in the laboratory frame, and in the 
limit /3 — >• (no precession) becomes the angle of rotation around the z axis. The 
magnetic moments canonically conjugate to the angles a, /3, and $ are P = niz — 
m^, mi3, and m^ respectively, where rriz is the z component of the magnetization 
m in the laboratory coordinate frame, m^ is the projection of m on the ^ axis 
of the rotating coordinate frame (see figure^, and mp is the projection of m on 
the line of nodes N, which is perpendicular to the axes z and ^. The free energy 
density consists of three terms, F = Fq + F\r + V, where 



m 



Fo = m H 

2x 



(48) 



includes the magnetization and the Zeeman energies, and the gradient energy Fy 
and the order-parameter dependent energy V (dipole energy in the case of ^He) 
will be determined later on. Here x is the magnetic susceptibility. 

For the phenomena observed experimentally only one degree of freedom is essen- 
tial, which is connected with precession, i.e., with the conjugate pair "precession 
phase a-precession moment P" . In contrast to the mode connected to the longitu- 
dinal magnetic resonance (oscillations of the longitudinal spin component), which 
was discussed in previous sections, the precession mode is connected with the trans- 
verse magnetic resonance (nuclear magnetic resonance in the case of ^He), in which 
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ruz does not oscillate essentially. 

Neglecting nutation, the directions of the axis ^ and of the moment m coincide. 
Then /3 is constant, m^ = m, mp = 0, P = m{cos f3 — 1), and the free energy 
density Fq becomes 

2 
77? 

Fo = — -{m + P)H, (49) 

where H is a strong constant magnetic field parallel to the z axis. The Hamilton 
equations for the precession mode are: 

da 5F d(Fs/ + V) 

9^=^5P=-"^ + ^ dP ' 
IdP 6F 0(Fv + F) 9Fv 

7 ot da oa oVia 

where col = jH is the Larmor frequency and F = f d^RF is the total free energy. 
Since this section deals with nuclear spins, here 7 is the nuclear gyromagnetic ratio. 
In the experiment the magnetization amplitude is determined by the magnetic field: 
m = xH- 



Equation (50) describes free precession without dissipation. One of the most 
important mechanisms of dissipation is the Leggett-Takagi mechanism [22j related 
with the process of equilibration of the magnetization of the normal component 
with the precessing magnetization of the superfluid component (for details see the 
reviews by Bunkov [13] and Fomin |12j and references therein). This mechanism 
becomes ineffective at low temperatures, and this leads to the Suhl instability of 
the uniform precession, which is discussed below. So one cannot observe uniform 
precession in ^He-B at very low temperatures. The dissipation leads to a precession 
decay, and in order to support the state of uniform precession in the experiment 
the energy dissipation must be compensated by the energy pumped by the rotating 
transverse magnetic field. Assuming that the balance of the pumped energy and 
the dissipated energy eventually leads to stationary precession, we may further 
ignore the both. The stationary precession state corresponds to the extremum of 
the Gibbs thermodynamic potential, which is obtained from the free energy with 
the Legendre transformation: G = F + ojpP/^, where the precession frequency 
Up = —da/dt plays the role of the "chemical potential" conjugate to the precession 
moment density P. 

Up to now the theory was rather general and valid for precession in any magnet- 
ically ordered system. Referring to ^He-B particularly, the dipole energy in ^He-B 
is [Ellis] 



V 2X^ 



2 r in2 

-— ^ (1 + cos<I>)u + cos$- - 
157^ [ 2 



(51) 



where is the longitudinal NMR frequency and u = cos/3. At the stationary 
precession the angle $ does not vary in time and can be found by minimization of 
the Gibbs potential. In the state of uniform precession without spatial gradients 
only the dipole energy depends on $, and the equation for <I> is 



dV Ax^ 



2 



d^ 1572 



'1 + cos <^)u + cos $ - - 



2 



(l + u)sin$ = 0. (52) 
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Solution of this equation yields 



cos$ = ^'y^ — -, y(u) = Ofor /3< 104° (u> -1/4), 
8x^^ /I 



cos$ = l, V{u) = -^[-+u] for /3 > 104° (u < -1/4). (53) 



Thus at n > -1/4 (/3 < 104°) and at u < -1/4 (/3 > 104°) one must choose two 



different branches of the solution of equation (52). The critical angle j3c = 104° is 
called the Leggett magic angle. 

It has already been known about 50 years from studies of nonlinear ferromagnetic 
|36j and antiferromagnetic [37] resonance that the state of uniform spin precession 
with finite precession angle can be unstable with respect to excitation of spin 
waves (Suhl parametric instability). Though Suhl instability is a phenomenon of 
the nonlinear classical wave theory [38] it is easier to qualitatively explain it in 
terms of spin- wave quanta (magnons). The processes leading to instability are 
transformations of n quanta of uniform precession into two spin- wave quanta with 
wave vectors itfc: 

rujjp = u}{k) + (jj{—k). (54) 

The precession is unstable if at least one of these processes is allowed by the laws 
of energy and momentum conservation. The three-magnon process (n = 1) cor- 
responds to the first order Suhl instability. The process is possible if a quantum 
of uniform precession of frequency ujp can dissociate into two quanta of the lower 
spectral branch with frequency u)p/2. Another possibility to destabilize uniform 
precession is a four-magnon process of transformation of two quanta of uniform 
precession into two quanta of the same spectral branch with finite wave vectors itfc 
(the second-order Suhl instability, n = 2). The process becomes possible if a non- 
linear correction to the frequency of uniform precession has an opposite sign with 
respect to the frequency dispersion du"^ /dk'^. In the theory of nonlinear waves the 
latter condition is called Lighthill's condition, which is necessary for modulation 
instability [39]. The second-order Suhl instability is an example of it. In all known 
examples of uniform precessions in magnetically ordered systems there are condi- 
tions for at least one type of Suhl instability. In superfluid ^He Suhl instability is 
possible in the A [401 [4lJ and B |32t H3] phases. But as any parametric instabil- 
ity, the Suhl instability can be suppressed by dissipation, which leads to a critical 
precession angle below which the state of uniform precession remains stable. In 
the B phase stable uniform precession is possible at temperatures T > 0.4 Tc. At 
lower temperatures dissipation is weak and cannot block the Suhl instability. This 
explains a sudden transition to the regime of "catastrophic relaxation" observed 
by Bunkovet al. g^Q 



7.2. Stability of spin-precession supercurrents (Landau criterion) 

Let us consider the state with uniform spin-precession current proportional to the 
gradient Va. The total free energy should now include the gradient energy of ^He-B 



^In contrast to Surovtsev and Fomin j42i i43j . who explained catastrophic relaxation with the bulk Suhl 
instability, Bunkov et al. 145] suggested another mechanism of Suhl instability, which exists near the surface 
(see arguments over the two mechanisms by Fomin ;46 and Bunkov et al. |47 p 
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2 7^ 



,V$^ 



+ ci 



Vu' 



2(1 



(55) 



where 



X 



A{u) = A,[4{1 - uY + ciil - u% 



T 



(56) 



n = cos/3, and cy and c± are longitudinal and transverse spin wave velocities. The 
expression for F^ assumes that all gradients are normal to the axis z parallel to 
the dc magnetic field. 

An important feature of the dipole energy in ^He-B is its independence of the 
precession angle a. In accordance with Noether's theorem this means that the 
precession moment is strictly conserved. Thus the equations describing stationary 
precession with frequency cop are 



iOp 



-UJL +7 



5{Fy + V) 
6P 



(57) 



V- J = 0, 



(58) 



where 



dF 

dVa 



-A{u)'Va 



(59) 



is the spin-precession current. 
Apart from Va, other gradients are absent: Vu = V$ = 0. Then equations (53) 



and (|55|)-(|57|) yield the following equation for u: 

^2 



(cop - col)col - [4 + (ci - cJ)'u]Va^ = at /3 < 104° [u > -1/4), 



(cop - col)ool - [cl + (ci - 4\)u]Va^ + ^ ^~^ '^' = at ^ > 104° {u < -1/4). 



15 



(60) 



The proper way to check stability of the current state given by equation (60) is to 



use the Landau criterion |35j . Since we check stability of the relative minimum at 
fixed averaged gradient Va, we must do a new Legendre transformation choosing 
a new Gibbs thermodynamic potential 



G = G + J ■'Va = F + cop— + J • Va, 

7 



(61) 



which has a minimum at the specified values of the precession P = M(l — uq) and 
the gradient Vao. Now we must find the energy increase due to fiuctuations. It is 
easy to check that fluctuations of $ always increase G, so it suffices to retain in 
the fluctuation energy only terms quadratic in small deviations u' = u — uq and 
Va' = Va — Vao from the stationary values uq and Vao (we omit hereafter the 
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5G = A{u 



(Va' 



,n2 



2 +A'{u)u'Va'+^ 



V"{u)+A"{u) 



Va' 



,'2 



(62) 



The fluctuation energy is positive definite, i.e., the current state is stable, so long 
as Va does not exceed the critical value 



VOr 



A{u)V"{u) 



1/2 



4f] / cj(l-n)2 + ci(l-n2) V/' 



A'^{u) - A"{u)A{u)/2\ 5^/3 I [cj + (c^ - c2)n]2 + c\/3 



(63) 



This expression yields the critical gradient on the order of the inverse dipole length 
1/^^ = 0/c_L, but it is vaUd only at n < -1/4. At u > -1/4 {(3 < 104°) the dipole 
energy vanishes, Voc, and the superfluid precession transport is impossible. 

Another definition of the critical gradient was suggested by Fomin |39]: He be- 
lieved that the spin-current state can be stable as far as the gradient does not 
exceed the value Voc = -s/iiop — ujl)ujl/c_\_, which is the maximum gradient at 



which equation (60) for u has a solution. Fomin's theory allows stable supercur- 
rents for /3 < 104° when the dipole energy and the Landau critical gradient vanish. 
Arguing in favor of his critical gradient, Fomin f50] stated that the Landau crite- 
rion is not necessary for the superfluid spin transport since emission of spin waves, 
which comes into play after exceeding the Landau critical gradient, is not essen- 
tial in the experimental conditions (see also the similar conclusion after equation 
(2.39) in the review by Bunkov [l3]). This argument is conceptually inconsistent. 
If the experimentalists observed "dissipationless" spin transport simply because 
dissipation were weak, they would deal with ballistic rather than superfluid trans- 
port. An example of ballistic spin transport will be considered in section [8} As was 
stressed in section [l} the essence of the phenomenon of superfluidity is not the ab- 
sence of sources of dissipation, but ineffectiveness of these sources due to energetic 
and topological reasons. The Landau criterion is an absolutely necessary condition 
for superfluidity. Fortunately for the superfluidity scenario in the ^He-B, Fomin's 
estimation of the role of dissipation by spin-wave emission triggered by violation 
of the Landau criterion is not conclusive. He found that this dissipation is weak 
compared to dissipation by spin diffusion. But this is an argument in favor of im- 
portance rather than unimportance of the Landau criterion. Indeed, spin-diffusion, 
whatever high the diffusion coefficient could be, is ineffective in the subcritical 
regime, in which the gradient of the "chemical potential" is absent. On the other 
hand, in the supercritical regime the "chemical potential" is not constant anymore 
and this triggers the strong spin-diffusion mechanism of dissipation. 



7.3. Experimental evidence of the superfluid spin-precession transport 

As was discussed above, the appearance of spin current itself is not yet a manifes- 
tation of spin superfluid transport. Supercurrents appear in spin waves or domain 
walls where they transport spin on distances of the order of wavelength or width 
of domain walls, but hardly it would be reasonable to call it superfluid transport. 
Similarly spin currents in the domain wall separating HPD from the bulk without 
precession cannot be a manifestation of superfluid transport. A convincing evidence 
of spin superfluidity would be spin transport on long distances. This evidence was 
presented by Borovik-Romanov et al. [51j studying spin current through a long 
channel connecting two cells filled by HPD. The schematic set up of their exper- 
iment is shown in figure l9l There is a dc magnetic field parallel to the vertical 
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Figure 9. Spin-precession transport through a channel connecting two cells 
filled by HPD. The horizontal arrow shows the direction of the spin-precession 
current in the channel. The precession angle /3 in the channel is less than in 
the cells (the analogue of the the Bernoulli effect, see the text). 



axis z. The HPDs in the two cells were supported with independently rotating rf 
magnetic fields and with different precession phases as a result of it. A small dif- 
ference in the frequencies of the two rf fields leads to a linear growth of difference 
of the precession phases a in the cells. This creates a phase gradient Va in the 
channel accompanied by a spin-precession supercurrent. The rf coils can monitor 
precession phases in different parts of the set up. Due to a linear growth of Va 
in time eventually it reaches the critical value at which a 2tt phase slip occurs. It 
is possible to register this event via its effect on NMR absorption [T3]. Thus the 
critical gradient can be measured as a function of the precession frequency. 

Despite aforementioned conceptual fiaw of Fomin's theory, Borovik-Romanov et 
al. |5HI52] found an agreement of this theory with the experiment. Let us compare 
now the experiment with the theory based on the Landau criterion. Equation (63) 
for the Landau critical argument contains the value u = cos/3 inside the channel, 
which is different from u in the cells where there is no precession phase gradients 
Va (see figure ^. This is an analogy with the Bernoulli law in hydrodynamics 
(liquid density is less in areas with higher currents) . The value of u in the channel 



grows with Va according to equation (60) at fixed precession frequency ujp. The 



latter is controlled in the experiment and in stationary states does not vary in 
space, exactly like the chemical potential in stationary states of superfiuids. If 



(wp - U}l)u}L < 



A'(-l/4)7%^„ ,2 



X 



-(Va, 



(64) 



where Va, is the critical gradient from equation (63) at u = —1/4, u reaches the 



value -1/4 earlier than Va reaches the Landau critical gradient. Since no stable 
supercurrent is possible at u > —1/4 the critical argument is determined as a 
solution of equation (60) at n = —1/4: 



Va, 



5c2 - c2 



1/2 



(65) 



The experiment was done at small ujp — ul, and its results must be compared with 
equation (65). For the ratio Cn/c^ = 4/3 [50j the latter gives the value of Va, 

by the numerical factor ■y/12/17 ~ 0.84 smaller than Fomin's result, which was in 
about 1.5 times larger than the critical gradient in the experiment. Thus the theory 
based on the Landau criterion even better agrees with the experiment |5H [52] , and 
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an approximate agreement of Fomin's result with the experiment cannot be used 



as an argument in its favor. At larger values of ul — ojp, when the condition (64) 



is not satisfied, the critical argument is determined by (63) and not proportional 
to \/{ojp — ujl)^l- So the difference with Fomin's theory becomes more essential 
(for more details see reference [53j). 

Further important development in experimental studies of spin-precession super- 
fluidity was observation of a spin-current analog of the Josephson effect |54j . The 
weak link was formed by making a constriction of the channel (orifice) connecting 
the two cells. 



7.4. Spin-precession vortex and its nucleation 

As was already discussed in section [4j at gradients less than the Landau critical 
gradient, the barrier, which impedes the current decay, is related to vortex motion 
across the flow streamlines (phase slips). Similarly, vortices called spin-precession 
vortices appear in the spin-precession flow, and the vortex core radius was estimated 
to be on the order of the dipole length ^^ |35j. The barrier for vortex growth in 
the phase-slip process vanishes at phase gradients of the order of the inverse core 
radius. So the threshold for vortex instability agrees with the critical gradient from 



the Landau criterion [equation (63)]. This is usual in the conventional superfluidity 
theory [TT] . 

Later Fomin |50j showed that the vortex core must be determined by another 
scale ^F = c±/ yj{ujp — ujl)^l, where up and ujl are the precession and the Larmor 
frequencies. This was supported by Misirpashaev and Volovik f55] on the basis of 



the topological analysis. According to equation (60) in the ground state without 
spin currents (wp — u}l)^l = 16r2^|ti + 1/4|/15. So if u is not too close to -1/4 
S^d = c_l/0 and ^p are of the same order of magnitude. But if n — t- —1/4, i.e. the 
precession angle /? approaches to the critical value /3c = 1.82 rad (or 104°) the core 
radius becomes Vc ^ ^p ^ Cd/{P — f^c), i-e., by the large factor l/(/3 — /3c) differs 
from the earlier estimation re ~ ^d |35]. So the latter is valid only far from the 
critical angle, where /3 — /3c ~ 1- Since no barrier impedes vortex expansion across 
a channel if the gradient is on the order of 1/rc, the large core re ~ ^d/{(3 — Pc) at 
(3 ^- (3c leads to the strange (from the point of view of the conventional superfluidity 
theory) conclusion: The instability with respect to vortex expansion occurs at the 
phase gradients ~ 1/rc essentially less than the Landau critical gradient ~ l/Cd, 
obtained for any (3 > f3c. Recently a resolution of this paradox was suggested [56] : 
At precession angles close to 104° at phase gradients less than the Landau critical 
gradient but larger than the inverse core radius, no barrier impedes phase slips at 
the stage of vortex motion across streamlines, but there is a barrier, which blocks 
phase slips on the very early stage of nucleation of the vortex core. So for these 
gradients stability of current states is determined not by vortices but by vortex-core 
nuclei. 

It should be stressed that, in contrast to the previous subsection, where the 
growth of Va was accompanied by the growth of u at fixed {ujp — u}l)ul, the 
present analysis is performed at /? = arccosn fixed in the channel excepting an 
area a vortex core or its nucleus. It never exactly equal to f3c though (3 — /3c could 
be whatever small. Thus (ujp—ujl)ul grows with Va. Vortex nucleation starts from 
a " protonucleus" , which is a slight localized depression of the superfluid density 
[determined by A{u) in our case]. The nucleus, which is related with a peak of a 
barrier, corresponds to an extremum (saddle point) of the Gibbs potential given 



by equation (61). Therefore, the nucleus structure should be found from solution 



of the Euler-Lagrange equations for this Gibbs potential. The first step is to vary 
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the Gibbs potential with respect to a. Let us restrict ourselves with a ID problem, 
when the distribution in the nucleus depends only on one coordinate x. Then the 
distribution of Va is given by Va = —J/A{u), where J is equal to the spin- 
precession current J = —A{uoo)'Vao, which is determined by the gradient Voq far 
from the nucleus center. Expanding with respect to small deviation g = u — Uoo 
from the fixed value Mqo at infinity one obtains 



G 



J2 



2A{u) 



Xc\ 



i^g)' 



72 2(1 
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J2 
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+ V{u, 
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+ V{u) 
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+ 



g^ J2 (fA{i 



6 2 



dv? 



(66) 



where we took into account that d^V {u) / dv? = 0. The linear in g terms must 
vanish at the stationary current state. The term quadratic in g determines the 
stability of the current state: it vanishes at the Landau critical current 



J. 



d''V{u^) \d^[A{u^y 



dv? 



dv? 



(67) 



This is exactly the Landau critical argument Vqc, which was determined in sec- 
tion 7.2 [equation (63)]. Considering the case of the current close to the critical 
value and using the Taylor expansion of A~^{u) around u = —1/4 one obtains 



G 



1572 



2 2 



9 



(68) 



where 



a = 0.239 



2 \ 2 

xciJ 



j2), 6 = 0.577 



\xcl) 



J 



The Euler-Lagrange equation for this Gibbs potential, — A^r = ag 
determines the distribution of g: 



(69) 
bgV2 = 0, 



9 = 50 1 - tanh 



(70) 



where go = 3a/b = 1.24(j2 - 

2/^ = A.lXclh'y^ 



r, 



j'^)/J^ is the value of g in the nucleus center and 
- J2 is the nucleus size. The energy of the nucleus. 



1572 



3a/b 



a9 



bg^ 



dg 



64xc2 



,5/2 



2572 62 



-5 = 0.2145 



( J2 _ j2)5/2 



J'c 



(71) 



determines the barrier for the process of the vortex core nucleation. Here S is 
the cross-section area of the channel. Since in the limit J ^ Jc the nucleus size 
rp diverges, our ID description is always valid close enough to the critical point, 
where r^ ^ \/S. When r^ becomes smaller than the transverse size of the channel, 
one should consider the 3D or 2D (in the case of a thin layer) nucleus. The first 
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stage of this problem is to find the distribution of Va from the continuity equation 
V[5fVa] = 0. Its solution demonstrates that outside the nucleus the distribution 
of Va is the same as around the vortex ring (3D case) or the vortex dipole (2D 
case). In particular, in the 2D case 



Va = Vao 



g{ri)rldri 



Vao 2r(r • Vqq) 



(72) 



In contrast to the ID case, the relation between Va and g is not local, so the 
following variation of the Gibbs potential with respect to g leads to an integro- 
differential equation. However using the scaling arguments one may conclude that 



the nucleus energy can be roughly estimated from the expression ( 71 ) for the ID 



case with replacing S by r^ or by VpL for the 3D and the 2D cases respectively (L 
is the thickness of the 2D layer) . 

This analysis demonstrates an unusual feature of the superfluid spin-precession 
transport at the precession angle close to the critical angle 104°: A bottleneck of 
the phase slip process is not connected with expansion of already formed (i.e., with 
sizes exceeding the core size) vortices but with the early stage of vortex nucleation. 
The spin-precession vortex in ^He-B was detected experimentally [571158). 



8. Ballistic spin transport by magnons 

Another interesting case of spin transport in magnetically ordered system is con- 
nected with magnons [591 [60] . This is an analogue of the normal mass current in 
superfluids, which arises due to transport of mass by quasiparticles (e.g., phonons 
at low temperatures). Let us find first what contribution to the spin current comes 
from one magnon. The simplest case is a magnon in an isotropic ferromagnet. 
Metastable spin-current states are impossible in such a ferromagnet but it is not 
essential for now: We look for a "normal" spin current. 

We consider an isotropic ferromagnet subject to a magnetic field H with the free 
energy density 



F = -M ■ H + -ViM ■ ViM. 



(73) 



The free energy is invariant with respect to rotations around the magnetic field H 
(the z axis), so in accordance with Noether's theorem the z component of spin is 
conserved. 
In the homogeneous ground state the spontaneous magnetization Mq is par- 



allel to H. Linearizing the Landau-Lifshitz equation (14) with respect to small 
perturbation m = M — Mq (m X Mq) one obtains 



1 dm 



-aVi[Mo X Vjm] + [H xm] 



(74) 



This equation yields spin waves ex e* "^* with the spectrum 



(75) 



which in contrast to the spectrum ( 19 ) without magnetic field, has a gap ^H equal 



to the frequency of the ferromagnetic resonance. 
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In the linear approximation the spin current J^ , which is determined by equa- 



tion (18), vanishes after averaging over the period of the wave. So we need the 



terms of the second order in m, which yield the current 

J^ = a{z-[mxVjTn\) =akj{m^). (76) 

The energy density in a single spin- wave mode is e = {uj/2^MQ){m'^)V, where V 
is the sample volume. Contributions of a magnon with the energy e = hco to the 
squared transverse magnetization and to the spin current are 5{m?') = AfiBM^/V 
and (5J^ = hv{k)/V respectively. Here hb = 7^/2 is the Bohr magneton and 
v{k) = duj/dk is the magnon group velocity. If there is an ensemble of magnons 
with the distribution function n{k) the total spin current will be 



Y,n{k)v{k) = ^ j dkn{k)v{k). (77) 



For the axisymmetric equilibrium Planck distribution the spin currents vanishes. 
But the spin current appears if the magnon distribution is the Planck distribution 
with a drift velocity Vn'- 

™o(^' ^) = eftKfc.g)-fc..„]/T_i - C^^) 

The magnon drift velocity i)„ is an analogue of the normal velocity in a superfluid 



liquid. Expanding the right-hand side of equation (77) in small Vn one obtains the 
linear relation between the spin current and the drift velocity. The Plank distribu- 
tion with a drift is valid only if interaction between magnons is more effective than 
interaction of magnons with other quasiparticles (electron, phonons), or lattice de- 
fects. Otherwise the magnon distribution must be determined from Boltzmann's 
equation, and the spin current appears only if there is a gradient of the magnetic 
field H, which plays a role of the chemical potential for spin since H = —dF/dm. 
The spin current proportional to the gradient of H is accompanied by dissipation 
and is determined by "spin-conductivity" J^ /VH. 

Another way to obtain the spin current in the magnon system is to connect a 
quasi-one-dimensional ferromagnetic channel of finite length with two bulk ferro- 
magnets (they act as reservoirs for spin) via ideal contacts [60]. Magnons cross 
the channel without scattering (ballistic regime). A spin current appears if there 
is a difference /S.H of the magnetic fields in the two leads. The physical picture 
is similar to that for electron ballistic transport analyzed in the framework of the 
Landauer-Biittiker approach [61J: The "right-movers" (magnons moving to right) 
and "left-movers" (magnons moving to left) are described by the equilibrium dis- 
tribution inside the leads, which magnons come from. Bearing in mind that the 
wave vector has the only component along the channel the spin current is given by 



r 
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' f°° dio 




/ dk „ 


'I'K 


[Jo dk 



/O j^ 

dk—no{k,H) 
-oo dk 



(79) 



Here no{k,H) is the Planck distribution given by equation (78) at i>„ = 0, taking 
into account the dependence of co on H. This yields the spin conductance J^ /AH 
obtained by Meier and Loss |60j . The origin of dissipation is similar to that for 



ballistic charge transport along a ID channel [6T]: There is no dissipation in the 
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channel itself, but inagnons arriving to the leads have the distribution different 
from the equilibrium distribution inside the lead. Its relaxation to the equilibrium 
leads to dissipation. 



It is worthwhile to stress again (see also section 7.2) that though dissipation in 
the channel is absent in the ballistic regime, this is not superfluidity: In the ballistic 
regime dissipation is absent because there is no sources of dissipation. 



9. Equilibrium spin currents in helimagnets 

Though the previous analysis addressed the case of ferromagnets it can be ex- 
tended on anti- and ferrimagnets. The analysis is also relevant for spin transport 
in spinor Bose condensates of cold atoms, for which the stability of spin-current 
states (spiral structures) was also analyzed in the spirit of the Landau criterion 
[52] . The condition for superfluid spin transport is a proper topology of the mag- 
netic order parameter: the magnetic order parameter space can be mapped onto 
a circumference with only weakly broken symmetry (or no broken symmetry at 
all) with respect to rotation around the circumference. However, spin currents in 
helimagnets still require a special discussion. 

The magnetic structure in helimagnets is a spatial rotation of the magnetization 
(spin) in the easy plane. This structure appears due to Dzyaloshinskii-Moria in- 
teraction linear in gradients of magnetization, which break invariance with respect 
to space inversion. Its energy is given by D ■ [M x [V x M]] |17] . This energy 
is relativistically small and can affect only the direction of the magnetization M, 
but not its absolute value M. In easy-plane magnets M is fully determined by its 
angle in the easy-plane, and the free energy is 

F= j d^RF= I d^R I ^ + ^^^^^ + K[l - cos(ny)] - DM'^V .fX , (80) 

where it is supposed that the magnetic spiral is oriented along the axis z. This 



energy differs from the free energy in equation (33) with the term linear on the 
angle gradient V z'f. 

The term linear in the angle gradient does not affect the equations of motion 
but it is important for definition of the equilibrium magnetic structure: the phase 
gradient is present in the ground state. It also changes the expression for the spin 



current replacing equation (26) with 



J^ = -^ = -A(V^-A.), (81) 

where k = (DM'^/A)z. The ground state is determined by minimization of the free 
energy with respect to the average gradient (Vip), which determines the density 
n{V(p)/2Tr of domain walls. Focusing on the limit of low density of domain wall 



see equation (41)] the free energy density is 



^(Vv^) ( -— =^ - A; ) . (82) 



Thus at A; < 4:/^:^/nl (strong anisotropy) there is a complete "phase fixation", 
spin being directed along some of the in-plane easy axes, and the ground state is 
uniform (V(/9 = 0). This means that the spin current J^ = Ak is present in the 
ground state. At A; = 4/Try/n there is the phase transition to the spiral structure 
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with (Vif) 7^ 0, which is a chain of sohtons of the sine Gordon equation obtained by 
variation of the free energy. Near the phase transition the equihbrium spin current 
is stiU close to J^ = Ak. However in the hniit k ^ 1// the anisotropy energy is 
not essential and V(y9 ~ k. Then the equilibrium spin current is absent despite the 
presence of the spiral structure with Vc/3 7^ 0. 

The equilibrium spin currents in helimagnets were calculated by Heurich et al. 
|63j and by Bostrem et al. [M]. In their calculations rotational invariance in the easy 
plane was broken by an inplane magnetic field H±. This corresponds to the free 



energy (80) with K = MqH± and n = 1. Our analysis shows that an equilibrium 
spin current is a generic feature of any helimagnet even without an in-plane field 
since one cannot imagine an easy-plane magnetic material without at least some 
finite in-plane anisotropy. The phase transition between the phase-fixed state and 
the the state with a helical structure is a typical example of the commensurate- 
incommensurate phase transition, which is common in condensed matter physics 
|65j . In particular, such a transition driven by a in-plane magnetic field was studied 
experimentally [66j and theoretically [67J in the quantum Hall bilayers. 

In addition to equilibrium spin currents there are also possible metastable spin 
currents. Their stability is determined like it was done in section |4j but in the 
Landau criterion, Vc^ < -s/Ea/A, and in the expression, which determines the 



barrier for vortex expansion [equation (32)], one should replace the phase gradient 
V99 with |Vc/? — k\. 

The helical structure at the equilibrium (and possibly in the metastable state) is 
stationary and does not move. However, as discussed in section |6j spin injection to 
a medium, which does not support superfluid spin transport, is impossible without 
some nonequilibrium magnetization. Then the helical structure must move, but this 
motion can be blocked by pinning. The effect of pinning on moving incommensurate 
structures is well known [65]. It is similar to the effect of coercivity on motion 
of domain walls. However, in the limit of high spin currents most effective for 
dissipationless spin transport pinning is suppressed due to strong overlapping of 
domain walls. If they overlap so strongly that the phase (angle) gradient does not 
vary in space, the pinning force must vanish. 

Equilibrium spin currents are possible not only in helimagnets. They were re- 
vealed also in spin-Heisenberg rings subject to inhomogeneous magnetic fields 
|68t I69| . These currents are connected with the geometric phase quantization in 
ring geometry and, being inversely proportional to the length of the ring, are pos- 
sible only in mesoscopic rings. Their origin and properties are similar to mesoscopic 
currents arising in systems without magnetic order, which will be considered in sec- 
tion [H^i 

Equilibrium spin currents were considered by a number of researchers as a bizarre 
and unphysical outcome. They believed that they are background currents, which 
do not correspond to real spin transport |70H72j . This stimulated attempts to 
redefine the definition of spin current in order to avoid spin currents in the ther- 
modynamic equilibrium. For Heisenberg magnets it was done by Schiitz et al. |72j . 
In many aspects confusion about a proper spin-current definition in magnetically 
order systems reflects the similar problem in systems without magnetic order, and 
we postpone its discussion for section [TT} 

10. Magnon Bose-Einstein condensation vs spin superfluidity 

Connection of Bose-Einstein condensation (BEC) with superfluidity was debated 
from the very beginning of the theory of superfluidity. This connection is not 
straightforward. BEC not necessarily leads to superfluidity as the simple case of 
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the ideal Bose-gas demonstrates. On the other hand, one may explain superfluidity 
without referring to the concept of BEG as Landau did creating his two- fluid theory. 
There are debates also what the term BEC really means and how widely one may 
use it. Nowadays there is a strong tendency to see manifestation of BEC in various 
cases of formation of coherent states in magnetically ordered systems. The BEC 
scenario was suggested for the phase transition to the antiferromagnetic state from 
the saturated ferromagnetic state in a strong magnetic field in magnetic insulators 
TlCuClOa [73l [71], and CS2CUCI4 [75]. In yttrium-iron-garnet films Demokritov 
et al. [76] [see also [77]] observed condensation of magnons in a single state in- 
terpreting it in terms of BEC. Finally, Bunkov and Volovik |78H80] ) reinterpreted 
spin-precession superfluidity in ^He-B (see section u\ in terms of BEC. Dealing 
with semantics it is difficult to decree strict rules and to reach a broad consensus. 
Without trying to do it, we shall simply follow recent debates on the question when 
one may and when one may not use the label "magnon BEC". 

I start from a general argument why the term "magnon BEC" is at least problem- 
atic (but not forbidden!). Let us remind what did BEC mean at the time of Bose, 
Einstein, and London when the concept was coined. The quantum theory showed 
that if identical particles are described by the Bose-statistics, at low temperatures 
they can condense at the state of the lowest energy. Since all Bose-condensed par- 
ticles are described by the same wave-function, the BEC yields a new coherent 
classical field, which cannot be imagined in the classical theory of particles. So 
from this point of view the non-trivial essence of the BEC is appearance of a co- 
herent classical field, which does not exist in the classical description of particles. 
Returning now to the magnon BEC, a magnon can be a result of quantization of 
a classical field (spin waves). If one claims magnon BEC, this in fact means that 
"particles" magnons can be described by a classical fields. However, this is a logical 
circle: First we state that we need to quantize a classic field and to introduce quanta 
of this field, then we tell that we should describe condensed "particles" (quanta) by 
a classical field. So we return exactly to the same classical field, which we started 
from! In other words, in this special case quantization of spin waves was not neces- 
sary. Calling coherent spin waves by a Bose-condensate of magnons, one must allow 
also to call laser modes "Bose-condensates of photons" , or deformed states of solids 
"Bose-condensates of phonons" . Analogy between Bose-condensates and coherent 
states was well known long ago |81j . But analogy does not mean identity. Some 
reasonable restriction should be imposed on using the term BEC. It seems that 
more careful is to use this term only referring to objects, which can be particles 
in the classical limit. This means that coherent sound (phonons), electromagnetic 
(photons) or spin (magnons) waves are not Bose-condensates, but excitons in semi- 
conductors are. However, no clever classification can avoid "gray areas", where it 
fails to provide definite conclusions. Can one call coherent polariton states as Bose- 
condensates keeping in mind that the polariton is a superposition of an exciton and 
a photon? I leave this question without an answer. 

Putting aside this argument as too rigorist let us discuss various claims of 
"magnon BEC". They can be classified as equilibrium phenomena, or non- 
equilibrium ones, which require energy pumping for their observation. Let us start 
from the equilibrium case of the phase transition to the antiferromagnetic state 
from the saturated ferromagnetic state in a strong magnetic field in magnetic insu- 
lators TlCuClOa tZaiZl], BaCuSi206i82j, and CsgCuCU [75j. The theoretical basis 
for it was known from 80s |83j (see also [H]): The transition from the ferromagnet 
(all spin parallel to the magnetic field) to the two-sublattice antiferromagnet with 
spins of sublattices in the plane can be described in terms of magnon Bose-operators 
in the ferromagnetic state. This transition is one from numerous examples of ori- 
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entational phase transitions in magnetism. If the BEC label is used for this case 
it could be used for any other. More generally, any transition described in terms 
of soft-mode instability (e.g., structure transitions in solids) may be considered as 
BEC. Then the phenomenon of BEC loses its uniqueness. Reservations concerning 
using the term "magnon BEC" for equilibrium phenomena were expressed in a 
number of publications [73 ES] • 

Another class of magnon BEC is non-equilibrium states, which require en- 
ergy pumping for their support. There are two ways to support a coherent non- 
equilibrium state. The first one is to pump energy at the same frequency with 
which the phase of the coherent state varies in time. An example of such non- 
equilibrium coherent states is the coherent precession in ^He-B discussed above in 
section [7j In these cases there is a magnetic field rotating with the same frequency 
as the precession angular velocity. Snoke [86] called these examples of coherent 
states "driven condensates" suggesting that in such cases the word "condensate" 
was not appropriate and contradicted to the spirit of the term "condensate" as a 
spontaneous thermodynamic phenomenon. Putting aside this argument for a while 
and accepting the view of Bunkov and Volovik [78H80J . who considered the co- 
herent spin precession as a BEC phenomenon, the question arises why they do it 
only with respect to the coherent precession in ^He-B. The coherent precession 
takes place in the nonlinear ferromagnetic or antiferromagnetic resonances, which 
were studied more than 50 years. Already at that time it was well realized that 
precessing magnetization may be treated purely classically as a gigantic magnetic 
dipole, which emits coherent radiation similar to Dicke superradiance [87j|. Later 
experimental evidence of this radiation was reported [88]. Following the semantic 
criteria of Bunkov and Volovik one may call it "discovery of magnon BEC" . 

Really spontaneous magnon condensation was observed by Demokritov et al. [76] 
in yttrium-iron-garnet films: they pumped microwave energy at higher magnon 
modes and observed accumulation of magnons in the lowest-energy magnon state, 
which corresponded to nonzero wave vector. It is important that pumped magnons 
were incoherent ^89j, so ensuing accumulation in a single state was spontaneous in- 
deed. Earlier the phenomenon was discussed theoretically by Kalafati and Safonov 
[SD]. Condensation requires some threshold pumping level. Accepting undisputed 
similarity of this phenomenon with BEC, it is worthwhile to draw attention to 
another analogy, which is not less "intriguing" than the BEC analogy: the analogy 
with lasers. Keldysh |91| defined "lasing" in his review on the exciton BEC (see 
p. 255 there): "The accumulation of a macroscopic number of initially incoher- 
ent excitation quanta in a single-photon mode is lasing" . Replacing "single-photon 
mode" by "single-magnon mode" one can apply this definition to the experiment 
of Demokritov et al. |76| . On the basis of this analogy and since magnons are cou- 
pled with electromagnetic waves, one should expect coherent microwave radiation 
similar to laser radiation. This radiation was really observed recently by Dzyapko 
et al. [92]. Radiation results from confluence of two magnons with opposite wave 
vectors (one- magnon radiation is forbidden by the momentum conservation law) . It 
was demonstrated that the frequency of radiation is determined by the condensed 
magnon frequency but not the frequency of the pumping signal. Description of 
magnon accumulation in a single state in terms of coherent magnon states was 
suggested many years ago [93] by analogy with coherent photon states in lasers. 
On the basis of this idea recently Rezende |94j developed the theory of coherent 
microwave radiation detected by Dzyapko et al. [92]. If it were possible to check 
coherence of this radiation, it would be a demonstration of magnon laser. 

In order to distinguish between BEC and laser, it was suggested to use BEC 
trademark only if there is a quasi-equilibrium distribution of non-condensed 
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magnons with a well defined chemical potential that requires an effective magnon 
thermalization [79j (in lasers non-condensed photons are not in equilibrium) . As one 
of possible semantic restrictions on using the term BEC, this criterion of BEC is for- 
mally legitimate. Demokritov et al. pS] experimentally demonstrated that magnon 
condensation in the momentum space in yttrium-iron-garnet films was accompanied 
with the quasi-equilibrium Bose-Einstein distribution of magnons. On the other 
hand, since the property of coherence is related only with condensed magnons it is 
difficult to understand why the distribution function of non-condensed magnons is 
so crucial. 

Up to now we discussed the property of coherence but not of superfluidity. It is 
necessary to stress that coherence does not lead automatically to superfluidity as 
sometimes assumed [80j. An example of the coherent state without superfluidity 
was presented in section Hi an (anti)ferromagnet fully isotropic, or with uniaxial 
easy-axis anisotropy. There is long-range correlation of spins in these cases, but 
no metastable spin-current states. The property of superfluidity is determined by 
topology of the order parameter for a concrete type of a condensed state. Let us 
check it for various BEC cases considered in the present section. The equilibrium 
phase transition in strong magnetic fields from ferromagnetic to antiferromagnetic 
ordering certainly leads to potential possibility of superfluid spin transport: In 
the antiferromagnetic phase the strong magnetic field keeps the magnetizations 
of the sublattices in the plane normal to the field. So spiral structures with the 
antiferromagnetic vector rotating in the plane can be stable and provide dissipa- 
tionless spin currents. The precession states in ^He-B also can be accompanied by 
stable spin-precession currents as was demonstrated theoretically and experimen- 
tally (section It]). As for magnon condensation in yttrium- iron- garnet films [76], one 
cannot expect superfluid spin transport in this case since there is no easy plane 
for order-parameter rotation in the experimental geometry (the magnetic field was 
parallel to the film). However, if the modification of the experimental geometry (the 
magnetic field normal to the film) recently suggested by Tupitsyn et al. [96] were 
realized, the question on possible superfluid spin transport would become relevant. 



Part II: Spin currents without magnetic order 



11. Definition of spin current 

Up to now we considered spin currents, which are possible only in magnetically 
ordered media. Meanwhile, a lot of attention was devoted to other types of dissi- 
pationless spin currents, which can appear without magnetic ordering. They are 
analogues of persistent charge currents in normal metals [97j and can appear even 
in equilibrium. Certainly in this case they cannot relax and are genuinely persistent, 
in contrast to superfluid spin currents, which are only metastable. 

We shall analyze the spin transport in normal systems using the simplest model: 
noninteracting electrons. Then one can start from a single electron in magnetic 
(H) and electric (E) fields, which are not uniform in general. The electron wave 
function is a two-component spinor 



*=('$;')■ (83) 



March 3, 2010 17:4 Advances in Physics SpinRev 



Advances in Physics 37 

We need the Hamiltonian with relativistic corrections |98] : 

where <I>o = hc/e is the single-electron magnetic- flux quantum and cr is the vector of 
Pauli matrices. In vacuum the constant A is equal to e/4mc^, but in crystals it can 
be much larger [14j . The electric field E may include not only an external field but 
also crystal fields. The crystal electric field is possible in 3D systems with broken 
space-inversion symmetry. In 2D electron systems on semiconductor surfaces the 
effective electric field normal to the surface can result from asymmetry of the 
electrostatic potential confining the 2D system to the surface (structure inversion 



asymmetry). This leads to Rashba spin-orbit coupling (section 13.1). Since the 
Pauli matrices do not commute, the spin-orbit term proportional to E can be 
considered as a non-Abelian gauge field similar to the Yang-Mills fields j99j . 



Writing down the Schrodinger equation for the Hamiltonian (84), 



5* 

ih— = n^, (85) 



one can derive the continuity equations for charge: 

O(^tijr) 



dt 
and for spin: 

dSp _ ;ia(*t(j^^) 
~dr ~ 2 dt 

with the following expressions for the currents: 



V • j, (86) 



V-/ + G^, (87) 



m 1 2^ ^ $n 



/ = ^ |-^(*V^V* - V*t^^) _ ?^(*t^^) 
2m I 2 cPo 

--[{^^a/sa + aa^)^) x E]\ . 



(89) 



Since spin is not conserved, there are source terms (torques) in the spin-balance 
equations: 

G^ = - Y[^ X (*ta*)]^ _ ^{*t[^ x[Ex{V- 27rM/cI>o)]]/3* 

+ [[E X (V + 27riA/$o)*"^] x ct]^^}. (90) 



The spin current, equation (89), agrees with the current definition used by Rashba 
|100j and many others: 

j^ = -{^^{apv + vap}^), (91) 
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where 



v = — [ -iV - — Act X E] 92 

m\ $0 7 

is the operator of the electron group velocity. 

One may write down the spinor wave function ^ as 

* = ^(sin|e^^O' ^''^ 

where n = (^t^) is the electron density, 7 is the tipping angle of the average spin 
with respect to the z-axis, 0-|- and (pi are the phases of the two spinor components. 
Note that here we consider a single-electron wave function, and the electron density 
is normalized to unity: / n{R) dR = 1. Introducing the global phase cj) = \{(t)^ + (j)i) 
and the relative phase (pz = <P\.~ '^t ^^^ averaged spin depends on the latter, which 
is the angle of rotation around the z axis in the spin space: 

, . (*ta,*) . , / , (*Vj/*) . . , 

(o-x) = = sin7cos02, {ay) = = sin7sin0z, 

n n 

{'yz) = - ^^ = cos7. (94) 

n 

In the variables n, 7, (p, and (j)z the electron (free) energy is given by 



F= / SRF 



. 2 7 (^^ 2ttAV 




2I (^. 2^^V 
cos — Vfflt z — 

2 V ^ $0 y 



fisB cos 7 H [z X EjVc, 

2m 



f^2y2 ^ 

- [(cr) X £] • ( V0t + V04) + ^^E^ + ^^ (Vn, V7) \ 



(95) 



where the two last terms are not essential for the further analysis. The charge and 
the z spin component currents are: 

e dF neh f_ , 27rA cos7_ , ,„ , ^,1 ,^^. 



' C0S7 V(/) - — — ^ -\[z>iE]\ = -\^3 + J^ 



dV4)z 2m \ ' y -r ^^ j ^ l jj 2e 

(97) 

Here we divided the spin current on two parts. The first one is a convection: the 
charge current is transferring also the average spin. The second one may be called 
a pure spin current, which is determined in the coordinate frame moving together 
with electrons, i.e., with the velocity v = j /en: 

J' = l^n {-^V<A. - m - {a){az)) ^ E]\ . (98) 

Zm, Z 
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Neglecting the term proportional to the electric field this spin current is propor- 
tional to the gradient of the angle 0^ of the spin rotation around the axis z in 
accordance with Noether's theorem. The coefficient before the gradient (stiffness) 
is proportional to the squared inplane component of the spin (oc sin'^ 7) similarly 
to the spin current in magnetically ordered media (see the first paragraph of sec- 
tion |4]). 

Let us discuss various terms in the expressions for charge and spin currents. The 
terms related to the global-phase gradient V(/> and the electromagnetic vector A are 
common and do not require any special comment. The terms proportional to the 
spin-angle gradient ^(pz are related with the geometrical (Berry) phase |101H103| . 
which results from the transport of spin over a closed path in a parametric space, 
which is the configurational space in our case. For spin ^ the Berry phase is —2^) 
where SI is the solid angle subtended by the varying spin during the cyclic process. 
So the Berry-phase "gradient" is VF = — | Vfi = — 1(1 — cos j)^(j)z- The gradient 
is in quotation marks because it is a well defined gradient of a scalar function only 
at constant tipping angle 7. In general the Berry phase depends on the path used 
for its definition. In order to separate the Berry-phase current from that connected 
with the global phase the latter must be redefined. Introducing (j) = cp^ instead of 
(/) = 5((/>t + 4>i) the charge current is 

j = -n\v^-Vr--^-X[{a)xE]]. (99) 

It is worthwhile to stress that the Berry phase not only path-dependent but is not 
single- valued even for a given path. Indeed, one may redefine the solid angle ri in the 
Berry phase F = — 2^ replacing dQ = 27r{l — cos j)d(j)z by dO = — 27r(l + cos7)(i</>j;. 
So the Berry phase depends on what pole of the sphere was included into the solid 
angle. The difference between two values can be accounted for in the phase </>. Thus 
the Berry phase is a part of the global geometric phase, which the electron obtains 
after moving around the closed trajectory, and separation of the Berry phase from 
the rest part of the global phase is not unique. 

The charge and the spin currents also contain a spin-orbit term proportional to 
an electric field. These terms are a manifestation of the Aharonov-Casher effect 
|1Q4] related with the non-relativistic interaction of the magnetic moment (spin in 



our case) with the electric field. The Rashba spin-orbit interaction (section 13) and 



the spin Hall effect (section 14 ) originate from this term. It can also be attributed 
to the Berry phase for the parallel transport in the momentum space [103] . 

We have discussed spin currents in the free-electron model, but the presence 
of various components of the current is determined by symmetry and topology 
without being restricted with a particular microscopic model. For example, the 
Aharonov-Casher effect was considered in the tight-binding model [105] and in the 
presence of random spin-orbit interaction |106j . As already stressed in the end of 
section [3l presumption that spin transport requires mobile carriers of spin [T31 [21] 
is not justified: For phenomenology it does not matter whether spin is transported 
by itinerant carriers or via exchange interaction between localized electrons. 

There were worries in the literature about ambiguity of the spin-current defini- 



tion [To]. Indeed, the definition of the spin current given above [equation (91)] is 
not the only possible choice. One can redefine the spin current by adding to it any 
current (5jf (jf — >• j^ +Sj^), if it is accompanied by redefinition of the spin torque: 
Gfs — )• G/3 + ViSj^ . This is a purely formal ambiguity of current definition (like 
freedom to choose various definitions of potentials in electrodynamics), which must 
not lead to any ambiguity in physical predictions. This only means that any defi- 
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nition of current is not complete without accompanying definition of spin torque. 
In principle, it is possible even to avoid the spin-current term at all, treating the 
whole current divergence Vijf as a part of the spin torque. This would exclude 
the word spin current from the scientific lexicon. As discussed in section [6j one can 
avoid using such terms as flow or current, and describe the same phenomena using 
only concepts of deformation, spin stiffness, or torque |1 1 ^1631 FlOTj . However, it is 
not necessary: In many cases the "current language" provides a very transparent 
physical picture of processes in various spin systems. 

There were also attempts to redefine the spin-current so that to eliminate the 
torque term G^ from the spin continuity equation making it to look as a continuity 
equation for a conserved quantity. For example, Shi et al. [21j (see also the discus- 
sion of their work by Bray-Ali and Nussinov ||108j ) used the definition of a torque 
dipole term by Culcer et al. jl07j rewriting the torque density as a divergence of a 
torque dipole density P^: 

G/3 = -V • P^. (100) 

The torque dipole density was included into the redefined spin current 7"^ = j^ + 



pp. Then the spin continuity equation (87) becomes 

dS. 



dt 



+ V • r'' = 0. (101) 



However, there is a payoff for this formally legitimate operation: The spin current 
becomes non-local since it is determined by an integral over the torque distribution 
in the whole bulk. Eventually any choice of the spin-current definition must yield 
the same predictions for observations if the whole procedure is correct. But using 
nonlocally defined spin currents it is more difficult to study effects related to local 
processes, which do not conserve the total spin. We shall return back to this is- 
sue discussing the observation of spin currents in the equilibrium Rashba medium 



(section 13.4) 



12. Equilibrium spin currents in one-dimensional rings 

12.1. A ID ring in a homogeneous magnetic field 

The simplest example of the equilibrium spin and charge currents is persistent 
currents in a one-dimensional ring |10H I109J . Let us consider first a ring of radius 
i? in a constant magnetic field H = Hz normal to the xy plane of the ring, 
which provides the magnetic flux <l> = ttR^H through the ring. We neglect the 
Zeeman energy since its effect (splitting of the Fermi energies for spins up and 
down) becomes important for the magnetic fields ~ ^qu/R much higher than the 
fields ~ ^o/R^ when the Aharonov-Bohm phase shift oc $ becomes essential. Here 
n = N/2ttR is the electron density for A^ electrons in the ring. We consider zero 
temperature. Introducing the azimuthal angle (p for the position of the electron in 
the ring, the Hamiltonian is 

s2 / .a ■- = 



where p = $/$o- Taking into account the periodic boundary conditions, the elec- 
tron spinor components are eigenfunctions of the angular momentum: ^-f-,4,((/') oc 
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g««0^ where s is an integer. Let us consider electrons with the same spin (up or 
down). The total energy of the Fermi sea and the ground-state persistent currents 
depend on whether the number of electrons s even or odd. If the number 2p+l of 
electrons is odd their energy is: 



E = ^j;;;mE('-p)"- (lo^) 



-p 



This energy corresponds to the ground state only if —1/2 < p < 1/2. In other cases 
the electron Fermi see should be shifted. For example if p becomes larger than 1/2, 
one must sum s from —p +1 to p + 1. The particle current at —1/2 < p < 1/2 is 
given by 

e dE , . eh enh 

This function must be periodically extended on any p beyond the interval 
(—1/2, 1/2). This yields the periodic sawtooth dependence of the current on p with 
the period 1 . If the number 2p of electrons is even one should shift this dependence 
with the half-period 1/2. 

Taking into account the both directions of spin, the result depends on whether 
the electron number N is even or odd. If N is even the numbers of electrons with 
spins up and down are even or odd together. Then the charge current obtained for 
one value of spin is doubled, whereas the spin current, which is proportional to the 
difference of charge currents for two values of spin, vanishes. On the other hand, 
if N is odd the charge-current dependence for one spin direction is shifted with 
respect to another. This lead to the persistent spin current j^ = ^(j+ — j-)/2e 
with the period 1, where j± are charge currents for two values of spin. The spin 
current jumps between values ±h'^n/2mR every half-period. 



12.2. A ID ring in an inhomogeneous magnetic field: Berry-phase currents 

Let us consider now a ring in an inhomogeneous crown-shaped magnetic field dis- 
tribution. This problem was considered in a number of papers starting from Loss 
et al. |101| . There is a constant z-component Hcosx and a rotating transverse 
component: H^ = ff sin x cos i;^ and Hy = H sinxsiiKJ) (crown-shaped field distri- 
bution) . Here x is the tipping angle of the magnetic field on the ring with respect 



to the z axis (figure 10). The Hamiltonian for a single electron with a position 



characterized by the azimuthal angle (p is 

^ ^ 2 J?2 \~^d6 ~^] ~/^-B-f^{sinx[cos(/)(T^ + sin(/)(7y] +cosx'5-^}- (105) 

This Hamiltonian can be exactly diagonalized |1021 IllOj . The eigenstates are 
spinors 



l*+(0) 



COS ^e^'*^ 



sm 






, ^ ,,,, , COS 3^6*'-^ \ / -sin2e^'<^ . ,^^^, 
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Figure 10. A ID ring in a crown-shape magnetic field. The solid angle CI 
subtended by the spin parallel to H determines the Berry phase F = —nCl. 



with the eigenvalues of the energy 



e± = E±T 



AE^ 



+ {^bH sinxY (1 - cos 7) . 



(107) 



Here I is the quantum number for the orbital moment, which should include also 
the Aharonov-Bohm phase shift, and 7 is the tipping angles for the average spin 
with respect to the z-axis determined by 



A^ 



cos 7 



V^A£;2 + 4(^35sinx)2' 



(108) 
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Two states correspond to two opposite directions of spin with angles 7 and 7 + vr. 
The energies 

^+ = ^ — ^ ~ hbH cos X, 

h'^d + i)^ 

^- = o P2 + ^'BH cos X, (109) 

are obtained taking into account only the Zeeman energy in the longitudinal mag- 
netic field Hz = H cos x, and 

AE = E^-E+= \ T, ' + 2ijbH cos x- (HO) 

The signs + and - label the states with spin up and down in the limit of the 
longitudinal magnetic field {x — ?■ 0). In the eigenstates the average spin rotates 
in space around the z axis with the same speed as the magnetic field. Because of 
inhomogeneity of the magnetic field the states are not eigenstates of the operator 
of the orbital moment but are mixtures of the states with two different moments. 
However, they are eigenstates of the total moment (orbital moment + spin) with 
the quantum number / + ^ since cylindrical symmetry is not broken. 

The charge and spin currents are obtained by summation over all / and two 
possible spin states: 



eh 



2-KmR'^ 



^|;+l[l_cos7(/)]} + j;{/ + l-^[l-cos7(0]}|,(lll) 



h^ 



A-KmE? 



^|;cos7(0-^[l 



cos7( 



- j; |(/ + 1) cos7(0 + ^[1 - cos7(/)]} I . (112) 

The second terms ex (1 — cos 7) in the sums originate from the geometrical (Berry) 
phase jlOlj . The simplest case is the limit of strong magnetic field, when the spins 
in two states are parallel or antiparallel to the magnetic field, and 7 = x does not 
depend on I. Then the Berry contribution to the spin current j^ is proportional to 
the total electron density, while the Berry contribution to the charge current j is 
proportional to the differences of the densities in the two spin states, i.e., to the 
spin polarization. 

Similar effects exist in the presence of a spatially rotating electric field (external 
or crystal), which generates an effective magnetic field via spin-orbit interaction 
|1U3[ IllUj . The equilibrium spin currents can appear without electron mobility 
when no charge current is possible. An example of it is the persistent spin current 
in a spin-1/2 Heisenberg ring |68[ 169): Similarly to the case considered above, the 
crown-shape effective magnetic field makes the spin to subtend a solid angle, which 
leads to the Berry phase contribution to the spin current. 

All effects discussed in the present section are purely mesoscopic: The derived 
charge and spin currents decrease with increasing size of the ring. The case of 
macroscopic equilibrium spin currents, which are finite in the thermodynamic limit, 



March 3, 2010 17:4 Advances in Physics SpinRev 

44 E. B. Sonin 

is discussed in the next section. 



13. Equilibrium spin currents in the 2D electron gas with spin-orbit 
interaction 

13.1. Currents in a uniform 2D gas with spin-orbit interaction 

Spin-orbit interaction allows to govern the spin transport by electric field, which is 
expected to be useful for applications [TU] . This explains a great interest to systems 
with spin-orbit interaction. A classical example of such a system is a 2D electron 
gas with the Rashba spin-orbit term (let us call it Rashba medium). Rashba |100] 
revealed that in the Rashba medium spin currents appear even at equilibrium. 
However, he qualified them as "not real", which cannot lead to transport and 
accumulation of spin. Their presence in the ground state was considered as an 
inherent problem in the spin current concept [70, tl]. In the present section I shall 
address equilibrium spin currents in the 2D electron gas with spin-orbit coupling 
and analyze whether these currents have something to do with spin transport. 
The Hamiltonian for a 2D electron gas with spin-orbit interaction is 

2m I 

+i/3(r)[*t(^^V^.* - ayVy^) - (ci^V^^t _ ayVy^^)^)^ , (113) 

where a{r) and P{r) are parameters of the Rashba and Dresselhaus spin-orbit 
interaction respectively. In general they may depend on the 2D position vector r. 
For simplicity, we shall concentrate the further analysis on the Rashba interaction 
(/3 = 0) postponing its extension on Dresselhaus spin-orbit interaction for the end 
of this subsection. The Rashba term oc a is a particular case of the spin-orbit- 



interaction term in the Hamiltonian equation (84), assuming that XE = az. The 



Schrodinger equations for components of the spinor ^ are: 

h"^ f I 2 I ^i^i ■ ^'^i Ida i da 

m \ 2 dx dy 2 dx 2 dy 

.^ ; h'^ f l„o , 9V't . d'ijj^ Ida , i da , , ,_ ,, 



The spin current and the torque in the spin-continuity equation (87) are 



• ^2 fc2 

if = -'—{^^apVi^ - V,*t^^*) _ ^(*t{^^[^ X z]i + [ax z]^ap}^m^) 



^^ = "^ {(*^^['^ "" [^ "" ^]]^*^) " (^[[^ "" '^ "" ct]^^^}^)}- (116) 

The Greek super(sub)script /3 refers to three components x, y, z in the 3D spin 
space, whereas the Latin super(sub)script i is related with the two coordinates x, y 
in the 2D electron layer. 
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a) 



b) 





Figure 11. The ground state of the Rashba medium, a) The case km > a, 
the Fermi sea (shaded blue) fills the upper (+) and the lower (-) band, b) 
The case km < a, the Fermi sea fills only the lower band. 



In the uniform Rashba medium eigenstates are plane waves given by spinors 



1 



\/2 V" 



Akr 



(117) 



where k = =Fie*'^, ip is the angle between the wave vector k and the axis x (kx = 
kcosip, ky = ksinip), and the upper (lower) sign corresponds to the upper (lower) 
branch of the spectrum (band) with the energies (see figure 11) 



h'^ fk 



m 



lb ak 



h\kl 



a 



2m 



(118) 



The energy is parametrized by the wave number k^^ which is connected with ab- 
solute values of wave vectors in two bands as k = |/co =F ck|. The eigenstates are 
spin-polarized in the plane with spins 



h[k X z\ 
2 k 



(119) 



parallel or antiparallel to the effective spin-orbit magnetic field. There is no spin 
component normal to the plane {z axis). The group velocities in two bands are 
given by 



hk 2a r^ , 

\ \zx s\ 

m m 



hko k 

m k 



v(k) = — + —[zxs] = '-^^'-^. (120) 

mm m k 

Spin torque in the eigenstates is absent, but there are inplane spin currents: 



where £ij is a 2D antisymmetric tensor with components £xy = 1 and £yx = — 1. 
The spin current does not vary in space since there is no precession of spin oriented 
along the effective spin-orbit magnetic field. The latter is constant for a plane wave, 
and there is no torque on the spin violating its conservation. 
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Though any eigenstate is spin-polarized, after averaging over the equihbrium 
Fermi sea (we consider the T = case) the total spin vanishes. But the total spin 
currents do remain. The Fermi energy is ep = fi?{k'^ — a^)/2, where km is the 
maximum value of Uq. In the case k^, > a, when the both electron bands are filled 
(figure fllk) , the inplane spin currents are sums of contributions from two bands 

mi- 



Ji 



Airm ^'' 



+ a] kdk + 



k„i+a 



+ a] kdk 



3fe2 



Qirm 



eij{l22) 



At km < Oi only the lower band is filled (figure [Tip) , and 



h 



Anm 



k,„+a 



- +a ] kdk = Eij 

2 / 47rm 



1,3 



+ kma' 



(123) 



Here we presented the calculation of equilibrium spin currents in the 2D gas 
with Rashba spin-orbit interaction at zero temperature. Recently Bencheikh and 
Vignale |113j performed a more general calculation taking into account temperature 
effects and Dresselhaus spin-orbit coupling. In the presence of Dresselhaus spin- 



orbit coupling equation (122) is generalized to 



Gttttt- 



Q(a2-/32)e.^. + /3(a2-/32)(^^).^.] 



(124) 



The existence of equilibrium spin currents is related with broken space inversion 
symmetry and is a generic phenomenon in systems with spin-orbit interaction re- 
lated to the non-Abelian gauge invariance |114| . In order to demonstrate that an 
equilibrium spin current is able to transport spin, we should consider nonuniform 
media. 



13.2. Spin currents in a nonuniform Rashba medium 

Let us consider a slightly modulated Rashba medium with the Rashba parameter 
varying in space as [112J : 



a{r) = uq + ai cos{p ■ r). 



(125) 



The eigenstates found above must be corrected using the perturbation theory with 
respect to ai: ^ = ^o + ^- Here ^o is the spinor for a uniform medium with 



a = ao given by equation (117). The equations for the first order correction ^' are 
(in components): 



pAe^l - ao[kK*{k)+pK*{pM[ = ^^ 
-ao[kK{k) +pK{p)]tp\ + pAe^l = ^ 



where 



kK*{k) + 
kK*{k) + 



PK*{p) 
2 

pK*{p) 



e''''"^ cos{p ■ r), 
e**'"^ cos(p • r), 



(126) 



Ae 



h"^ fp 



m 



+ pk ^aok . 



(127) 
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The solution of linear equations for ^' should be used for derivation of all relevant 
physical quantities (densities, torques, and currents). The general expressions are 
rather cumbersome. Moreover, the perturbation theory fails in the limit p — )• 0. 
Therefore, we restrict ourselves with the limit p^ k, oq. The torques and currents 
for inplane spin components are (only linear in ai terms are kept) 



Gi±ik) 



2aiaoh pjk 
m p^ 



p2^2 



sin(p • r), 



(128) 



jWik) 



ai 



h^ 



2m 



-Pi- 



^isPs 



p^ 



+ £ij =F 



4e,, 



jp-k) 

,2^2 



21 



p 



cos(p-r). (129) 



The torque and the current for the z-component of spin are given by terms of higher 
order in 1/p and vanish after integration over the Fermi sea. The integration of 
the torque and the current for inplane spin over the Fermi sea yields for the case 
km > ao- 



Gi 



1 



4^2 / G,+ {k)dk+ ^^^ 



,h^ 



1 

47r^ 



Gi-{k)dk 



"^^^U^«o + ^)sin(p.r) 
nm p \ 3 ' 



(130) 



ai 



h'' 



8^2 



m 



^isPs 



■?i = i^ / fM^) ^^ + i^ / •?'-5-(^) ^^ 



n + 



STTEijao 
p2 



k^ao + 



a: 



3\ 1 



cos(p-r). (131) 



If km < Oio: 



Gi = ^ I Gi.{k)dk 



-^^i^^M + f )sm(p.r). 



(132) 



Sirm 



■I] 



Pj- 



eisPt 



p^ 



n 



S-KEijap 

p2 



kmOiQ + 



1 

4^ 

7,3 

3 



jUk)dk 



cos{p ■ r) 



(133) 



The first term in the spin current, which is proportional to electron density n, is 
divergence-free, whereas the divergence of the second term does not vanish and 
compensates the spin torque in the spin balance. Thus the second term is respon- 
sible for spin transport from areas, where spin is produced {Gi > 0) to areas where 
spin is absorbed {Gi < 0). One may consider this as a manifestation of spin trans- 
port even though it does not result in spin accumulation. Thus equilibrium spin 
currents can transport spin, and an attempt to distinguish equilibrium (persistent) 
currents from transport spin currents [115J hardly would be reasonable. 

It is worthwhile to note that the spin current in a modulated Rashba medium is 
linear in the spin-orbit coupling constant ao, whereas the dependence of the current 
on a in the uniform Rashba medium is cubic. Apparently the linear dependence is 
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electric current 



Figure 12. Spin-dependent reflection of electrons from an ideal impenetrable wall. The electron from the 
upper band {k\_) is reflected either as an electron from the same band (fcj. ), or as an electron from the 
lower band (fc^). 

a general property of nonuniform media. In particular, the same dependence was 
predicted by Sablikov et al. |116j . who considered equilibrium spin currents along 
the interface between the media with and without spin-orbit coupling. 



13.3. Interference and torque at edges of the Rashba medium 

One cannot understand the physical meaning of spin currents without a clear pic- 
ture of what is going on at borders of a system with a bulk spin current. Suppose 
that the Rashba medium occupies the semispace x < while at x > spin-orbit 
interaction is absent. The two semispaces have also different potentials, so the 
Hamiltonian is 



n 



2m 



+y(r)*t^, (134) 



where 



V{r) 



at X < 
U at X > 



a{r) 



a at X < 
at X > 



(135) 



The electron states near the interface between two regions with different spin-orbit 
constants have already been analyzed earlier \in\ I118J . The spin currents in a 
hybrid ring consisting of parts with and without Rashba spin-orbit coupling were 
also analyzed by Sun et al. |119j . 

At X < one should look for a superposition of plane waves: one incident wave, 
which is coming from x = — oo, and two reflected waves (figure 12). For high- 
energy electrons with kg > a and the incident electron in the upper band, the 
superposition is 



* 






1 



„ik^ 



+ ri 



1 



e *^+"^ + r2 



1 






(136) 
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where k± = =F^e*'^*, ip± = arctan{ky / k±x) are the aziniuthal angles of the wave 



vectors in the plane xy, and k±x = \ {ko ^ a)^ — fe| are the x components of the 
wave vectors corresponding to states of the same energy in the upper (+) and the 

lower (-) band. At x > the wave function is evanescent: ^ = I , I e^^vV e~^^^ ^ 

where kh = \/2mU /It? — k^ + a"^. 

The wave superposition should satisfy the boundary conditions, which include 
continuity of the both components of the spinor and jumps of first derivatives of 



these components |117l I118J related with derivatives of a in equation (114): 



(9* 

dx 



5* 

dx 



-iaay'^ . 



(137) 



The expressions for the reflection coefficients are rather cumbersome in general, 
and we restrict ourselves with the case of an infinite potential step at the x = 
(C/, fcfc -^ oo). Then the spinor wave function at x = vanishes, and expressions for 
the reflection coefficients become simple: 



ri 



K+ 



3«(V--</'+) + l' 



^2 



2e*'^- cos (p+ 

gi(ip_-(p+) _|_ I ' 



(138) 



The relation between the angles ip^ and (/?_ is determined from the condition that 
scattering does not change the component ky = (fco — a) sin c/?_|_ = (/cq + «) sin ip- . 
For low-energy electrons ko < a all three waves in the superposition belong to 
the two parts of the lower band, either to the right {k > a) or to the left (k < a) 
from the energy minimum (figure 11). If the incident wave corresponds to the state 
with k > a, the superposition is 



* 



Jkyy 



^ik+x _|_ ^^ 



-ik-f-X 



^2 



„ik-X 



(139) 



where k.± 



le 



iip± 



(p± = aician.{ky / k±x) , and k±x = a/(o i k^Y — ky. The positive 



sign before k- in the exponent of the second reflected wave was chosen because the 
negative group velocity of this wave. Since the electron transport is determined by 
the group velocity, the latter should be directed from the boundary into the bulk 
even though the wave vector is directed to the boundary. The reflection coefficient 
are 



ri 



K+ 



gi{ip+-ip-) _ I 



^2 



K+ 



2e-"^-cos(p+ 



(140) 



Whereas in plane-wave eigenstates of the Rashba Hamiltonian the spin has no 
z component, the interference between the waves in the superposition leads to 
partial spin polarization along the z axis. For ko > a the oscillating spin density 



March 3, 2010 



17:4 



Advances in Physics SpinRev 



50 



E. B. Sonin 



Sz = {h/2)'^^(Jz'^ (Friedel-like oscillation) is given by 



m = \{ri{ 



[e "2*^+ + ne-^^'^i^ + r2[e-^(^++^-) 



+l]e-i(fc+x+fc-.)^ + r^r2[e*('^+-^-) + l]e'^^ 



h{s\n (fj^ + sin c^_) cos ip 



+ c.c. 



,sm{2k+a,x) 

1 + C0S(V3+ - if-) 

■sin(A;+x.x + k^^x) — s\n{kj^,j;X — k-xx)]- 



fl41) 



Exchanging + and — one obtains the spin density s-z{k) for the incident electron 
from the lower band. Similar expressions can be derived for the low-energy case 
k^ < a, when all waves belong to the lower band. 

The expressions given above are valid only if ky < k^, or sin (/?_ < | /cq — a|/(fco + 
a). At k- > ky > k+ the reflection of the incident electron from the lower band 
to the upper one is forbidden by the conservation law. But the contribution of the 
upper band into the wave superposition is still present in the form of the evanescent 
mode. The wave superposition in this case is 
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where K- 
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(142) 



(143) 



The z spin density for this wave superposition contains not only the interference 
contributions but also the contribution from the evanescent component oc e^^: 



s_,(fc) = ^^°^ '^- [e^p^ + cos(2A:_^.x) - 26^^ cos(A;_^.x)] 
Kq sm (/3_ 

_^ cobif- _ (-;.^_^^)(,Qs2^_ir ■j^/2A;_^x) -2eP''sin(_^/c3;)]. (144) 

Kq sm c^_ 

This expression is valid independently of whether the electron energy is high (/cq > 
a) or low {ko < a). 

All contributions to the z spin density are odd with respect to the sign of ky and 
vanish in the equilibrium state. But in the presence of the voltage bias along the 
y axis the distribution function also has an odd component, and spin polarization 
becomes possible. This leads to the edge spin accumulation (polarization), which 
is important for investigation of the intrinsic spin Hall effect (section 14.3). 



The wave interference near the edge leads not only to z spin polarization but 
also to the spin torque. The existence of this torque is required by the spin balance 



(87): If there is no current in the vacuum and there is a bulk current normal 



to the boundary, the presence of an edge torque is inevitable and its total value 
(integral over the whole edge area) must be equal to the spin current from the bulk 
independently of particular properties of the edge (an ideally reflecting wall in our 
case). Moreover, the integral edge torque should compensate the bulk spin current 
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not only at the border with the vacuum but also at the interface between the 
Rashba medium and a medium, which does not allow dissipationless spin currents. 
Indeed, in the latter case spin diffusion accompanied by spin accumulation is the 
only mechanism of spin transport, but in equilibrium no dissipative process is 
possible and the spin current must vanish at the interface. 

But the type of the edge does influence the spatial distribution of the torque. 
We shall derive this distribution for a simpler case km ^ a (figure [lib) when 
all expressions can be expanded in /cq- Here km is the maximum value of /cq cor- 
responding to the Fermi level. In this limit the main contribution to the torque 
originates from interference of the incident wave with the second reflected wave in 



the superposition (139): 



Gy+ik) = -^^Re |e-^fe+-+i'=-- (i _ ^;^_) ^^j . (145) 

A similar contribution Gy- comes from the conjugate superposition, in which the 
incident plane wave oc e~ corresponds to the wave number k < a with the 
negative group velocity. The subsequent integration over the whole Fermi sea in 
the lower band yields 



Gy{x) 



■Km 



L-'*2 ( r, ) 1) <^ i '^rn-^ 



1/13 2 2A 2 



(146) 



where pFq(ai, ..., a^; 6i, ..., 6^; z) is the generalized hypergeometric function |120j . 
The total torque over the whole bulk J_ Gy{x)dx = fi^a^km/^T^fn exactly com- 



pensates the bulk spin current [see equation ( 123 ) in the limit km <C a] 

At large distances from the border the torques for single modes oscillate fast, so 
the asymptotic behavior of the torque can be analyzed using the steepest-descent 
method. This yields the asymptotic torque at x — >• — oo: 

^^ = -/l:4;^^Rw^^^"K^-i)- ^'''^ 

This Friedel-like oscillation may be suppressed by disorder or electron-electron 
interaction, which were neglected in our analysis. 

In summary, we have obtained the following picture of spin currents and torques 
in the restricted Rashba medium. There is no spin torque inside the medium far 
from medium edges, but there is a constant spin current there. On the other hand, 
interference of incident and reflected plane waves leads to spin torques of opposite 
signs (source and drain of spin) near the two edges. The role of the bulk equilibrium 
spin current is to transport spin from the spin source near one edge to the spin 
drain near the opposite edge. 



13.4. Experimental detection of equilibrium spin currents 

The central question for understanding the physical sense of the spin current is 
how is it possible, if possible at all, to detect the existence of spin currents ex- 
perimentally. This question is especially acute for equilibrium spin currents since 
they do not lead to any spin accumulation. However, spin current leads to electric 
polarization, which might be detected via electric fields produced by it. Indeed, 
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Figure 13. The cantilever with the integrated Rashba medium (thick solid 
line) on it. (a) A rigid substrate. The red arrow (above the Rashba medium) 
shows the direction of the spin current. The green arrow (inside the substrate) 
show the direction of the orbital-moment current. The currents result in me- 
chanical torques itr at the edges of the substrate (b) The substrate is now a 
flexible cantilever. The edge torque at its free end leads to its displacement 
h. 



the spin currents is a flow of magnetic moments related with spin. According 
to classical electrodynamics |121j a magnetic moment m moving with velocity 
V creates an electric dipole p = [i) x m,]/2c. The terms niiVj in this expression 
are proportional to the spin currents j*-: rfiiVj = 7J*-. So the dipole moment is 

Pi = ^ijklJk/'^c = eijkejl./2mc^. But this relation does not take into account crys- 
tal effect, which can strongly amplify spin-orbit interaction comparing with vacuum 
electrodynamics. A more reliable definition of spin-orbit dipole moment is using 
the thermodynamic relation p = d{H)/dE applied to the Hamiltonian ( |84[ ) [70]. 
This yields pi = EijkXjj.- For the Rashba medi um w ith t he sp in-orbit constant a 
the spin current j'*- is determined by equations (122) and (123), A = da/dEz, and 
the electric dipole is parallel to the z axis. 

The electric fields induced by stationary spin currents in conducting media were 
discussed by Hirsch |122j and Sun et al. |123) . These fields were also discussed 
for magnetically order systems [60l [69l 11241 I125J . where there are no itinerant 
carriers. This phenomenon is an effect inverse to the spin Hall effect (generation 
of spin current by an electric field), which will be discussed in the next section 



14 The inverse spin Hall effect was observed experimentally by Valenzuela and 



Tinkham |126j . Though the experiment was realized for spin-diffusion currents but 
not equilibrium spin current discussed here, there is no evident reason why the 
origin of the spin current would be essential for the existence of the effect. In any 
case, this suggests at least a Gedanken experiment ioi detection of equilibrium spin 
currents: as well as the charge currents in currents loops, which do not lead to any 
charge accumulation, are detected by magnetic-field measurement, one may detect 
equilibrium spin currents by electric-field measurement even in the absence of any 
spin accumulation. 

Recently it was suggested |127] to detect an equilibrium spin current in the 
Rashba medium by measuring a mechanical torque on a substrate at edges of the 
Rashba medium caused by the spin current. If the substrate is flexible, the edge 



torques should deform it (figure 13), and measurement of this deformation would 



provide a method to detect equilibrium spin currents experimentally. An appropri- 
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ate experimental technique for such a measurement is already known: a mechanical 
cantilever magnetometer with an integrated 2D electron system |128j . Earlier me- 
chanical detectors were suggested for detection of non-equilibrium diffusion spin 
currents [129]. A mechanical stress produced by spin currents in mesoscopic struc- 
tures with collinear magnetic order was also studied by Dugaev and Bruno |130| . 
They called it the magneto mechanical effect. 

Derivation of the mechanical torque produced by bulk spin currents is based 
on the conservation law for the total angular momentum (the spin + the orbital 
moment) in the system "2D electron gas + substrate". The continuity equation 



(87) for spin must be supplemented by the continuity equation for orbital moment: 



^ + S/^jP = -Gp. (148) 

Here Lp are /? components of the orbital-moment densities and J9 is the orbital- 
moment flux-tensor. The torque Gr in this equation is the same as in the con- 



tinuity equation (87) for spin, but appears with an opposite sign. This provides 
the conservation of the total angular momentum S + L. We consider currents of 
y components along the axis x, so f3 = y and 7 = x. At the equilibrium the time 
derivatives of momenta are absent. In section [13.31 we saw that inside the Rashba 
medium there is a spin current but no torque, while at the edge there is an edge 
spin torque, which compensates the bulk current. According to the total angular 
momentum conservation law this leads to an edge orbital torque and to a flux of 
the orbital moment with a sign opposite to that of the spin current, as shown in 



figure 13 Since the 2D electron gas has no y component of the orbital moment, 



the whole orbital torque must be applied to an edge of the substrate. Now if the 



substrate is a cantilever rigidly fixed at one end (figure 13), the mechanical torque 
T = Jx = —jx = J Gy{x)dx will deform the cantilever, and the displacement of its 
free end can be measured. 

The fact that the spin current must be accompanied by an opposite orbital- 
moment current was already noticed by Sheng and Chang [131] and Zhang and 
Yang LI32J. Transformation of spin to angular momentum at an edge of the Rashba 
medium was recently discussed by Teodorescu and Winkler |133j . Since the coun- 
terflow of the spin and the orbital moment does not lead to any flow of the total 
moment, Zhang and Yang [132] have concluded that the spin current is not observ- 
able and cannot induce electric fields discussed in the beginning of this subsection. 
However, this conclusion ignores the fact that the spin and the orbital moments 
have different gyromagnetic ratios. Therefore though the flow of the total mechan- 
ical moment really vanishes, the flow of the total magnetic moment does not. In 
fact, the orbital moment current, which compensates the spin current, is not an 
obstacle but an instrument for spin-current detection as is shown here. 

For numerical estimation of the effect we shall use the maximal value of the 
spin-orbit coupling ah'^ /m = 6 x 10~^ eV cm = 10"^*^ erg cm = 10~^^ J m 
quoted by Rashba [lllj for InAs based quantum wells. When calculating the spin 
torque, it was simpler to consider the case km ^ a. But the mechanical torque 



reaches its maximum at km > a [see equation (122)] when the spin current is 



jx = —a^h'^/6-Km ~ 10 ^ erg/cm=10 ^ J/m. In order to estimate the displace- 



ment h of the cantilever end (see figure 13), we use the cantilever parameters from 



reference [134j : the length I = 120 ;um and the spring constant k = F/h = 86 /^N/m, 
where F is the force on the cantilever end. Using the theory of elastic plates |135j . 
one obtains that the torque r = —jx produces the displacement h = 3t/2A;L 
This yields h ~ 0.45 ;um. Less optimistic estimations of the spin-orbit coupling 
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{ah?/m = 3 X 10~^ eV cm for InAs, or 1.4 x 10~^ eV cm for GaSb |136j ) predict 
an order or more smaller displacements, but certainly measurable with the modern 
micromechanical technique. The torque can be enhanced and tuned by an external 
magnetic field. 
In the presence of the external magnetic field one should add the Zeeman energy 



Hbct • H to the Hamilton equation (113). The spin current in a single-electron 



state is determined by the spin, which is parallel or antiparallel to the "effective" 
magnetic field H—a^o[kxz]/TT acting on the electron. Here <^o = hc/e is the single- 
electron flux quantum. Integrating the single-state spin current over the whole k 
space and assuming that the Zeeman energy is much larger than the spin-orbit 
energy, one obtains the bulk spin current 

where ep is the Fermi energy. In the interval —fisH < ep < fj-sH electrons flll only 



the lower band [the lower sign in equation (149)]. Then in terms of the electron 
density n = m{tF + usHj/lnh? 

If ep > IJ'bH, electrons fill the both bands, the contributions from two bands to 
the spin current cancel each other, and spin current vanishes in our approximation. 
The present analysis ignores the effect of the electromagnetic vector potential on 
the electron momentum, but this effect (which deserves a special analysis) is absent 
for the inplane magnetic field Hx- 

It is worthwhile to comment that ambiguity of spin-current definition, which 
was intensively discussed in the literature, has no impact on the effect considered 
here. As discussed in the end of section [TH one may redefine the spin current 
jj by adding to it an arbitrary term (jj — )• jj + 5jj) but at the same time it is 
necessary to compensate it by redefinition of the spin torque (Gj — )■ Gj + Vj5jj). 
If the balance of the orbital part of the angular momentum is also considered, the 
definitions of the orbital torque and flux must be compatible with those of the spin 
part, in order not to violate the conservation law of the total angular momentum. 
Eventually whatever definition was used any correct calculation must predict the 
same observable effect (displacement of the cantilever). After we defined the spin 



current by equation ( 115 ), we are not free anymore in the choice of the definition of 
the torque and the current of the orbital angular momentum: The fiux of the orbital 
angular momentum in the elastic cantilever should be defined as J* = eimnXmTj' 
where T" is the elastic stress tensor. This choice looks most natural since it defines 
the mutual torque between the spin and the orbital moment as a derivative of the 
spin-orbit energy with respect to the rotation angle. 
Let us look at an alternative definition of spin current [M| ■, which was discussed 

The "spin-conserving" current T^{x) = j^{x) + P"{x) = 



11 



in the end of section 

3^{x)-\- j^ Gis{x')dx', which includes the dipole torque term P^(x), is constant, and 
it is difficult to notice in this picture that there is a process of angular-momentum 
transfer between spin and orbital degrees of freedom. This is due to a nonlocal 
character of the current T^: it controls only the global balance of spin. Globally 
no change of spin occurs in the sample, spin being generated at one edge and ab- 
sorbed at another. Meanwhile, exactly Zoca/ torques, which do not violate the global 
spin balance, are responsible for the angular-momentum transfer to the orbital sub- 
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system, which leads to the mechanical deformation discussed in this section. 

14. Spin Hall effect 

14.1. Phenomenology of spin Hall effect 

If an electrical current flows through a conductor with spin-orbit coupling, this 
can give rise to a spin current flowing normally to the direction of the electrical 
current. This effect predicted by Dyakonov and Perel [137 ] was later called spin 
Hall effect |138j . In contrast to the usual Hall effect, the spin Hall effect originates 
from the effective magnetic field produced by spin-orbit interaction and does not 
require an external magnetic field for its existence. Experimental evidences of the 
effect in bulk semiconductors [139J, in a hole [140] and an electron [141] 2D gas 
have already been reported. Observation of the inverse spin Hall effect [126J also 
provides evidence for the spin Hall effect since the direct and the reverse effects 
are connected with the Onsager relations. 

There are two possible origins for the spin Hall effect. The first one, which was 
discussed by Dyakonov and Perel |137j . is due to spin-dependent scattering on 
impurities and is called extrinsic spin Hall effect. But broken space inversion sym- 
metry also can allow an intrinsic effect, which is not connected with impurities 
directly, though can be strongly affected by them (see below). 

In general spin currents in the spin Hall effect are not equilibrium currents con- 
sidered in the previous section. They are accompanied by dissipation and require 
an energy input for their existence, despite that dissipation is related not with the 
spin current itself but with the longitudinal charge current. What does unite them 
with the spin currents discussed through the present review, is common controver- 
sies about their definitions and the impact of the absence of spin conservation. We 
consider the 2D electron gas assuming the presence of spin-orbit interaction and a 
weak electric field parallel to the axis y. Independently of the microscopic origin 
of the spin Hall effect, the broken space inversion symmetry allows the current of 
the z spin component along the axis x given by 

f, = asuE, (151) 

where asH is the Hall spin conductivity. In a uniform system with constant asH this 
current is constant by definition. On the other hand, at the sample border the spin 
current must vanish. Thus one should look for a spin current of another nature or an 
edge torque, which would compensate the spin current of electric origin. According 
to Dyakonov and Perel [137] (see also the recent phenomenological analysis by 
Dyakonov |142) ) another current is a dissipative spin diffusion current, i.e. the 
total spin current is j| + j^, where 

Jh = -DsV.S,. (152) 

This current is not constant because of inevitable longitudinal spin relaxation de- 
termined by the time Ti. Thus the continuity equation for the z spin component 
is 

^ = -V.j!,-| = I).V15.-|. (153) 

For the stationary process with dSz/dt = and under the condition that the total 
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current vanishes at the sample border x = one obtains 

3h = -Jle--"'% 5, = -^e--/^% (154) 

where Lg = \/DsTi is the spin-diffusion length, which has already been introduced 
in section|6J Thus the spin Hall effect leads to accumulation of the z spin component 
at the sample edge (called also spin orientation). The degree of spin accumulation 
is governed by dissipation parameters Ti and Ds- A reader can find a detailed 
discussion of physical mechanisms responsible for dissipation in the review by Zutic 
et al. [10]. 



14.2. Extrinsic spin Hall effect 

The extrinsic spin Hall effect originates from the spin-orbit interaction related to 
the local electric fields induced by impurities or defects. A well accepted model for 
this interaction [1^ is described by the Hamiltonian term 

He,tr = Xa-[pxVViir)], (155) 

where p is the electron momentum and Vi{r) is the impurity axisymmetric poten- 
tial, which depends on the distance r from the impurity. The effect of the spin-orbit 
interaction on elastic scattering was known long time ago |144| . The interaction 
leads to asymmetry of scattering (skew scattering), and the differential cross section 
depends on spin of an incident electron: 

a±{e) = ao{9)[l±S{e)], (156) 

where o"o(^) is an even function of the scattering angle 9 and determines scattering 
of unpolarized electron beams, while S{9) is an odd function of 9 called Sherman 
function. The Sherman function determines the skew scattering. The upper and 
the lower signs correspond to spins +1/2 and —1/2. 

Various types of impurities and of 2D quantum wells were discussed in the litera- 
ture (see, e.g., |145j and references therein). Here we present the simple but general 
discussion in terms of effective cross sections similar to that by Engel et. al. [146J. A 
straightforward and physically transparent method to calculate a bulk spin current 
is solution of the Boltzmann equation. Studying the spin Hall effect they usually 
used the quantum Boltzmann equation, in which the distribution function was a 
matrix 2 x 2 in spin indices |15l I147H149J . However, for the extrinsic effect under 
consideration the z spin component is a good quantum number, which is affected 
neither by external electric field, nor by spin-orbit interaction determined by equa- 



tion (155). Indeed, the effective spin-orbit magnetic field is normal to the 2D layer 
plane and scattering cannot lead to any spin-fiop. This means that off-diagonal 
elements of the spin-density matrix vanish. Two diagonal elements correspond to 
two distribution functions f±{k) = fo{k) + f'±{k) for electrons with positive (+) 
and negative (-) z spins. Here fo{k) is the Fermi equilibrium distribution function, 
which does not depend on spin and direction of the electron wave vector k. The 
Boltzmann equation for the non-equilibrium distribution function /^ (k) generated 



by a weak inplane electric field E is 

^E^^ = n,VF |[/±(0) - f±i<P')]a±{9) d9, (157) 
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where vp is the Fermi velocity, rii is the impurity density, and (p and cp' = (j) + 9 are 
angles between k and E before and after a collision. The solution of this equation, 
as one may check by substitution, is 

Here e = h?k'^/2m is the electron energy, and the relaxation time 

(159) 



riiVFia'^ + 0-2 



is determined by the transport cross section, 

at= f i{e){i-cose)de, (leo) 

and the effective cross section related to scattering asymmetry, 

a,= f I{9)S{e)smede. (161) 

The parameter asjut is called transport skewness |146j . 

Knowing the non-equilibrium distribution function one can calculate the charge 
current parallel to E and the z spin current transverse to E by integration over the 
momentum space and summation over two spin values. Let us consider the zero 
temperature limit when dfo{e)de = 6{e — ep)- The charge current is given by the 
usual Drude formula, 

j = aE = -E^, (162) 

while the spin current jz = ctshE is determined by spin Hall conductivity |146j 

e repCTs h as .^„„. 

zvr n at 2e at 

Here a is the ohmic electron conductivity. 

It is believed that the extrinsic spin Hall effect was detected by Kato et al. |139] 
who observed spin accumulation on edges of n-GaAs layers. Engel et al. |146j have 
found a rough quantitative agreement of the experiment with the theory presented 
in this subsection, though the signs of the effects were opposite. This disagreement 
remains unresolved (see discussion in |146j ). 



14.3. Intrinsic spin Hall effect and edge spin accumulation 

The intrinsic spin Hall effect does not rely on spin-dependent effects in scattering 
but is related entirely to the uniform bulk spin-orbit interaction, which leads to 
spin-orbit-split band structure. The effect was proposed by Murakami et al. [150j in 
bulk p-type semiconductors using the effective Luttinger Hamiltonian for holes and 
by Sinova et al. |151j in 2D electron systems with the Rashba spin-orbit interaction. 
In the literature there were debates on the strength of the intrinsic spin Hall effect in 
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the 2D electron gas. Relating the spin current to the acceleration of electrons by the 
electric field in a pure material, Sinova et al. |151j concluded that the intrinsic spin 
Hall effect is characterized by the "universal" spin Hall conductivity asn = e/Svr. 
Originally it was believed that this result is insensitive to a small concentration of 
impurities, despite the fact that a steady state of a conductor in an electric field 
is impossible without collisions, which compensate the monotonous acceleration of 
electrons by the electric field. Later on it turned out that in the Rashba medium 
even rare collisions should not be ignored and eventually the intrinsic spin-Hall 
conductivity at zero frequency vanishes in an infinite system (see discussion and 
references in the review by Engel et al. [14j). But this conclusion is valid only for 
Rashba spin-orbit coupling linear in the electron wave vector k, and the intrinsic 
spin Hall effect for Rashba coupling nonlinear in k, or in the Luttinger model 
for spin-orbit coupled systems |152j is not ruled out. The intrinsic effect in the 
medium with nonlinear in k Rashba coupling follows from the theory using the 
quantum Boltzmann equation for the spin density matrix |1491 1153| . In contrast to 
the extrinsic spin Hall effect, the intrinsic spin Hall effect is not possible without 
off-diagonal matrix elements since the z spin current is proportional to them. This 
is because the bulk spin-orbit interaction keeps the spin inside the xy plane, and 
only correlation between eigenstates of the Rashba Hamiltonian (like interference 
between plane-wave states near edges) leads to z spin polarization and z spin 
currents. 

However, the experimental detection of the intrinsic spin Hall is rather problem- 
atic. The experimental evidences of the spin Hall effect reported in the literature 
were based on optical measurement of the spin accumulated on the sample edges 
presumably resulting from the spin Hall effect. Meanwhile, spin accumulation (po- 
larization) near boundaries might not be so simply related to spin currents in the 



bulk as repeatedly pointed out [T5 | I143t I149J . As was shown above (section 13.1), 
the bulk spin current is not necessarily accompanied by spin accumulation. More- 
over, spin accumulation at sample edges is possible even without bulk spin current. 
It was demonstrated for the ballistic spin Hall effect |154H156] . when the electron 
mean-free path exceeds the sample sizes. Edge spin accumulation without bulk 
spin currents takes place also in the standard collisional regime when the electron 
mean-free path is much shorter than the sizes of the sample but still longer than 
the distance where interference of incident and refiected waves responsible for edge 
accumulation takes place |157j . Let us discuss edge accumulation without bulk spin 
currents in more details. 

Since our goal is to analyze the case without bulk spin current, we may use 
the standard linear-in-momentum Rashba Hamiltonian. Moreover, we do not need 
to deal with the quantum Boltzmann equation and restrict ourselves with the 
classic Boltzmann equation for two scalar distribution functions corresponding to 
two diagonal elements of the spin density matrix. Indeed, according to references 
|148ill49] . for the linear-in-momentum Rashba interaction bulk spin currents vanish 
together with off-diagonal matrix elements, to which they are proportional. 



In section 13.3 it was shown that the eigenstates of the Rashba Hamiltonian are 
partially z spin polarized near the boundaries, though in the equilibrium there is 
no total spin polarization after integration over the Fermi sea. But in the spin Hall 
effect there is a charge current (along the y axis at our choice of the coordinate 



frame, see figure 12), and the non-equilibrium distribution function has a compo- 
nent odd with respect to the wave vector component ky. First we shall consider the 
ballistic regime when the voltage drop V occurs at the contacts, and there is no 
electric field inside the sample. In the narrow interval of energies ep + eV > e > ep 
around the Fermi surface only left-moving electrons with ky > are present. They 
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are responsible for the edge accumulation of the z spin. Bearing in mind that 
eV = de = {h'^/m)k±Fdk, the two band contributions to the total spin density 
Sz{x) = S-^^zix) + S-z{x) are determined by integrals over the Fermi circumferences 
of the two bands with the Fermi wave vectors k±p = \km =F o|: 

'meVk±F f^'"^ 
'^'^''^ = A^^Wk: i «±.(fc)d^±. (164) 

Here km is the value of ko at the Fermi circumferences. The spin densities s±z{k) 
for the eigenstates, which correspond to the wave vectors k of incident waves, were 



found in section 13.3, The asymptotic behavior of the spin density is determined 



by the evanescent-mode contribution [see equation (144)] and at x — )• — oo is given 

by 



meV I a 1 1 



Sn^n V km k^pk-F \x\-^ 
The total accumulated spin is given by |156j 

S. = f" ,^(,),, = ^ (ln|^ - 1^) . (166) 

J-oo S-K^ha \ \km-a\ km J 

For further comparison with the collisional regime it is convenient to connect the 
total spin not with the voltage V but with the electric current, 

J = li^ X hi:^ ^* ^- > « , (167) 

-nhkm I 1 at fcm < a 

where the 2D electron density is n = (A:^ + a^)/27r at km > a and n = akm/T^ at 
km < ot. Then 

mj f^^ km + a _2a\ ^ f ^ at fc^^ > a _ ^^gg^ 



87rena V \km — oi\ km J I km at fc^ < a 



In the limits of weak (a — t- 0) and strong (a — t- oo) spin-orbit interaction this yields 
5*2 = (r?2j'/247ren)(a^/A;^) and 5^ = —mj/ATren respectively. At km = a there is a 
logarithmic divergence, which can be cut either by the sample size or by nonlinear 
effects. 

Let us switch now to the collisional regime. As was explained above, since the 
bulk spin current is absent, one may use the standard Boltzmann equation for two 
scalar distribution functions f±{k) = /o(e) + f±{k) for the two bands, where /o(e) 
is the equilibrium Fermi distribution function, which depends only on energy e. The 
stationary solution of the Boltzmann equation for the non-equilibrium distribution 
function /' in a weak electric field along the y axis is |159| 

, erEdfoik) hkm . ^, . .^„„. 

Here the relaxation time r for elastic scattering on defects is determined by the 
transport cross section and in general is a function of A;. In principle r should 
differ for two bands. But the difference vanishes for weak spin-orbit coupling and 
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further will be neglected. The functions /4 determine the electric current equal to 
j = e^£'rA;^/27rm, for km > en and to j = e^ETakml'^'T^ni for km < a. The z spin 



densities for the two bands instead of (164) are given by 



s±z{x) = ^ I sm.ip±s±z[k)dip±. (170) 



Tedious but straightforward integrations similar to those for the ballistic regime 
yield the total edge spin: 



S, 



mj k'^ + a'^ \nn2 2\ i 2\Jakn 



327r2en /c^ 



3(fc„ — a ) arctan 



2\/afcm(3fe^ + 2kmOt + "io?) 2 o 2 



for the high-energy case km > a, and 

mj 



km- a 
+ 7T{km + 3a^ 



(171) 



S, 



IGir'^enakm 

^akm{^G? + ^k^ + 4afe 



?>(o? — kl,) arctan 



a — kr, 



a + kj^ 



"* + Anakn 



(172) 



for the low-energy case a > km- 

When a — )■ 00 the difference between the ballistic and collisional regime vanishes. 
On the other hand, in contrast to the ballistic regime, in the collisional regime the 
accumulated spin remains finite even in the limit of zero spin-orbit coupling q — t- 
being equal to 

This paradoxical result is explained by the divergence of the width ~ l/{k-x — k+x) 
of the spin accumulation area in this limit. In the ballistic regime this divergence 
is canceled after summation over the two bands. However, our analysis is valid 
only if all relevant scales including l/{k-x — fc+x) are less than the electron mean- 
free path. When this condition is violated the spin accumulation should go down. 



Figure 14 shows the reduced total accumulated spin Sz = AT:enSz/mj for the 
ballistic (curve 1) and the collisional (curve 2) regimes as functions of the density- 
dependent parameter a/^/ym. Edge spin accumulation without bulk spin currents is 
possible for other types of spin-orbit interaction. In particular, Bokes and Horvath 
[160] considered spin-orbit interaction related to the edge potential V{x) confining 
the electron gas: Hgo = chect • [k x ^V{x)]. They obtained the edge accumulation 
linear in the spin-orbit constant oe, which is insensitive to the details of the edge 
potential. This allows to expect that the assumption of the hard-wall potential 
made in our calculation for the Rashba spin-orbit interaction also is not so crucial 
for the final outcome of the calculation. 

Originally the spin Hall effect was defined as an effect related to bulk spin cur- 
rents. The present discussion shows that spin accumulation is not really a probe of 
the bulk spin current: the former can be absent in the presence of the bulk current 
and can appear in the absence of the latter. Therefore the question arises whether 
edge accumulation without bulk currents may be called the spin Hall effect. A 
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Figure 14. The plot of the reduced total spin Sz = 'inenSz/mj as functions of a/^/wn. 1 - the ballistic 
regime. 2 — the coUisional regime. 



choice of terminology usually is a matter of convention, taste, or tradition. Edge 
spin accumulation and spin currents require the same symmetry, and one may call 
the edge spin accumulation without bulk currents the edge spin Hall effect. 

In order to compare the edge and the bulk spin Hall effects, we scale the latter 
using the "universal" spin conductivity asH = j^ /E = e/Svr, though in reality this 
is far from being universal |14j . Here j^ is the bulk current of the z spin. Assuming 
that at the edge the bulk spin current is fully compensated with spin diffusion 
current (section 14.1), the total accumulated spin is jzTi = eETi/Sn, where Ti is 
the longitudinal spin relaxation time. So ratio of the edge to the bulk spin Hall 
effect is ~ r/Ti. 

For comparison with the spin Hall effect observed in the 2D hole gas |140^ I158J on 
may use r ~ lOh/Ep = 20 /kmVp, n = 2 x 10^^ cm~^, and the accumulation area 
width 10 nm given by Nomura et al. |158j . Then the total spin accumulated due to 
the edge spin Hall effect at a — )■ is about 70 % of the experimental value. So the 
interpretation of this experiment in the terms of the bulk spin currents probably 
must be reconsidered even if the spin-orbit interaction for these materials should 
be described by a model more general than used for the present analysis. 

There are other cases of electrically generated edge spin accumulation without 
bulk spin currents inside the sample. Adagideli and Bauer |161j discussed edge 
spin accumulation governed by spin diffusion, which occurred within the distance 
of the order of the spin-diffusion length Lg from the interface between the media 
with and without spin-orbit coupling. In contrast, interference spin accumulation 
discussed above, is not related to any dissipative process and occurs at the distance 
on the order spin-orbit length, which was assumed to be much shorter than the 
mean-free path. It is interesting that the diffusion governed accumulation provides 
the total accumulated spin ~ eEr [7T] of the same order of magnitude as the 
interference mechanism. Since Lg must essentially exceed the mean-free path, the 
interference mechanism provides much higher spin density but in a much narrower 
layer. Another example of edge spin accumulation without bulk spin currents will 



be discussed in section 14.5 addressing the quantum spin Hall effect 
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14.4. Spin currents in spin Hall insulators 



As was already mentioned, though transverse spin currents induced by the spin 
Hall effect are dissipationless themselves, the whole process is accompanied by 
dissipation in the longitudinal channel and therefore requires an energy source. 
Murakami et al. |162| (see also reference |163| ) proposed the totally dissipation- 
less spin-Hall effect in insulators with spin-orbit interaction. In a band insulator 
at zero temperature and the Fermi level within the gap between conductance and 
the valence bands a weak electric field cannot generate an electric current, but 
may result in a transverse spin current, i.e., the spin conductivity asH in the spin- 



Hall relation ( 151 ) is finite being determined by topological invariants of the band 
structure. Materials with this property are called spin Hall insulators. The dissi- 
pationless spin currents in insulators, when the Fermi level is inside the forbidden 
gap, were revealed by numerical calculations using realistic parameters of semi- 
conductor band structures |164j . The spin Hall effect is possible not only in band 
insulators. Meier and Loss [60] considered the spin Hall effect in a two-dimensional 
Heisenberg model consisting of localized spin, in contrast to itinerant electrons in 
band insulators. 

The spin Hall effect in insulators has a very important feature discerning it 
from the spin Hall effect in conductors. Since there is no charge current in the 
former case, there is no Joule heating and no energy input [163]. As a result, no 
dissipation process is possible, and the problem reduces to the equilibrium problem 
of an insulator in an electric field. The bulk spin current cannot be compensated 
by the dissipative mechanism of Dyakonov and Perel |137j , and one should look for 
a Hamiltonian mechanism of absorption (generation) of spin near sample edges. 
So there is a close analogy between spin currents in the spin Hall insulator and 



equilibrium currents in the Rashba medium analyzed in section 13, Indeed, the 
Rashba spin-orbit interaction is connected with an internal or external electric field 
normal to the electron-gas plane. The electric field is not able produce a current 
since electrons are confined in a 2D layer. So the 2D gas is an insulator in the z 
direction, but this does not rule out a Hall spin current along the layer. 

The equilibrium character of spin currents in spin Hall insulator impose a serious 
restriction on methods of spin-current detection. In particular, spin injection into 
a non-magnetic material without spin-orbit interaction is impossible, since in this 
material spin can be transported only by a dissipative (diffusion) current, which 
must be supported by energy pumping. On the other hand, it is possible to extract 
spin from a spin Hall insulator putting it into a contact with a magnetically ordered 
medium supporting dissipationless spin transport. 



14.5. Spin accumulation and the quantum spin Hall effect in topological 
insulators 

Recently great attention was attracted to remarkable properties of topological insu- 
lators, in which the quantum spin Hall effect was predicted [165^ I166J . The hallmark 
of a topological insulator is a forbidden gap originated from spin-orbit coupling and 
the presence of topologically stable helical edge states. Topological criteria and 
classification for these materials have already been carefully analyzed [167^169] . 
The outcome of this analysis most important for the goals of the present review 
is illustrated in figure [T5| At two edges of the sample parallel to the electric field 
(we address the 2D case) each spin is able to move only in one direction, which 
is opposite for two spin directions. The edge states cross the whole forbidden gap 
connecting the valence-band and the conduction-band bulk continua. The system 
is time-reversal invariant, and two states with opposite directions of the spin and 
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Figure 15. Edge states in a topological insulator. Wide blue arrows show spin direction (spin quantization 
axis is not necessarily in the plane as in the figure). At the upper edge spins up are rightmovers while 
spins down are leftmovers. At the lower edge directions of motion are opposite. The ballistic edge channels 
short-circuit the bulk insulator with infinite resistance. Therefore the whole sheet is globally a ballistic 
conductor without electric field inside. The external voltage bias V drops inside electrodes (see voltage 
distribution in the lower part of the figure) . 



the wave vectors form a Kramers degenerate pair. Because of helicity (one-way mo- 
tion) the edge states are robust against elastic backscattering and therefore may 
be treated as bahistic. This means that the ballistic edge channels short-circuit the 
bulk insulator with infinite resistance, and the whole sample is globally a ballistic 
conductor without electric field inside. 

Thus, whatever the bulk spin conductivity asH = j^ /E could be, bulk spin 
currents are absent simply because the bulk electric field is absent. Nevertheless, 
the voltage drop between two leads (see the voltage distribution along the sample 



in figure 15) leads to edge spin accumulation similar to that considered in section 



14.3 for bulk conductors. However, mechanisms of the edge accumulation in two 



case are different. Whereas in bulk conductors the spin accumulation arose from 
interference of electron plane waves reflected from the boundary, in topological 
insulators spin polarization appears because the leftmoving and rightmoving edge 
states with oppositely directed spins have different densities proportional to the 
voltage bias. The latter directly follows from the Landauer-Biittiker approach. The 
leftmovers occupy states below the Fermi level in the right electrode, whereas the 



rightmovers occupy states below the Fermi level in the left electrode (flgure 15). 



The difference between the two Fermi energies is eV ^ and the edge spin density 
(spin per unit length along the edge) is determined by the number of states in this 
energy interval: 



S, 



ese 
hvpi 



■V, 



(174) 



where vpe = de{ky)/hdky is the Fermi velocity of the edge mode with the spectrum 
£{ky) and Se is the effective "spin" of the edge states, which may be different from 
ib/i/2 in general. The parameters of edge modes were calculated numerically and 
analytically using simple but reliable models |1671 1168) . The word "spin" is in 
quotation marks since it is actually a combination of the electron spin (which is 
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zizh/2 of course) and the orbital moment of the electron in the band. Though it is 
dangerous to jump to general conclusions on which angular momentum eventually 
is relevant for various experiments, one may definitely expect that electromagnetic 
experiments (Kerr or Faraday effect, induced electric fields) address the magnetic 
moment in the edge mode. 

For estimation of the effective spin related to the magnetic moment one can 
consider the model usually used for the topological insulator |1671 I168J : the edge 
state crosses the forbidden gap separating the conduction and the valence bands, 
which originate from s-type {I = 0) and p-type (/ = 1) atomic orbitals. The p- 
type orbital corresponds to quantum numbers j = 3/2 for the total moment and 
rrij = ±1/2 for its projection on the quantization axis. The Lande factor, 

-, , JU + 1) - l{l + 1) + s{s + 1) 
^" = ' + MW) ' 

is equal to the electron g factor g^ = 2 for the s-type orbital and to the factor 
gv = 4/3 for the p-type orbital. The edge mode is a superposition of the two states 
with the weights te and t^ respectively {te + t^ = 1). Then the effective spin in 



equation (174) is 



Using this value of the effective spin the accumulated magnetization is obtained 



from expression (174) by multiplying with the electron gyromagnetic relation. 

It is interesting to compare the edge spin accumulation due to the quantum spin 
Hall effect with the accumulation due to the intrinsic spin Hall effect in conductors 



(section 14.3). Introducing the average electric field E = V/L, where L is the length 



of the sample, equation (174) at Se ~ /i yields Sz ~ eEL/vpe- Comparing this with 
Sz ~ eEr for the accumulation in conductors one sees that this transforms to the 
accumulated spin in topological insulators after replacing the relaxation time r by 
the time of flight L/vpe through ballistic edge channels. 

Soon after the theoretical prediction the topological insulators were experimen- 
tally detected in the HgTe quantum well jl70j by studying charge transport. It 
was demonstrated that at the quantum well thickness exceeding the critical value 
6.3 nm there was an interval of gate voltages where the conductance reaches the 
value 2e^//i independently of the sample width W (see figure 15). This is a clear 



evidence of the ballistic transport through edge states while the main bulk is not 
conducting. The topological insulators states were also detected in BiSb |171j . BiSe 
|172] . and BiTe [173J compounds by the methods of angle-resolved photoemission 
spectroscopy (APRES). In the literature the topological insulator state is called 
also the quantum spin Hall state. It is worthwhile of mentioning that despite im- 
pressive experimental demonstration of the quantum spin Hall state the quantum 
spin Hall effect itself is still wants its experimental confirmation: A "smoking gun" 
of edge spin accumulation have not yet been reported. 



15. Conclusions 

The present review focused on four types of dissipationless spin transport: (1) Su- 
perfluid transport, when the spin-current state is a metastable state (a local but 
not the absolute minimum in the parameter space). (2) Ballistic spin transport, 
when spin is transported without losses simply because sources of dissipation are 
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very weak. (3) Equilibrium spin currents, i.e., genuine persistent currents. (4) Spin 
currents in the spin Hall effect. The dissipationless spin transport was a matter 
of debates during decades, though sometimes they were to some extent semantic. 
Therefore it was important to analyze what physical phenomenon was hidden un- 
der this or that name remembering that any choice of terminology is inevitably 
subjective and is a matter of taste and convention. The various hurdles on the way 
of using the concept of spin current (absence of the spin-conservation law, ambi- 
guity of spin current definition, etc.) were analyzed. The final conclusion is that 
the spin-current concept can be developed in a fully consistent manner, though 
this is not an obligatory language of description: Spin currents are equivalent to 
deformations of the spin structure, and one may describe the spin transport also 
in terms of deformations and spin stiffness. 

The recent revival of interest to spin transport is motivated by emerging of spin- 
tronics and high expectations of new applications based on spin manipulation. This 
is far beyond the scope of the present review, but hopefully the review could jus- 
tify using of the spin-current language in numerous investigations of spin-dynamics 
problems, an important example of which is the spin Hall effect. 
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